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Cn Abstract. We undertake a comprehensive and rigorous analytic study of the 

evolution of radial profiles of covariant scalars in regular Lemaitre— Tolman-Bondi 
dust models. We consider specifically the phenomenon of "profile inversions" 
in which an initial clump profile of density, spatial curvature or the expansion 
scalar, might evolve into a void profile (and vice versa). Previous work in the 

^J\ literature on models with density void profiles and/or allowing for density profile 

inversions is given full generalization, with some erroneous results corrected. We 
prove rigorously that if an evolution without shell crossings is assumed, then only 
the 'clump to void' inversion can occur in density profiles, and only in hyperbolic 
models or regions with negative spatial curvature. The profiles of spatial curvature 
»' I follow similar patterns as those of the density, with 'clump to void' inversions only 

l^j\ possible for hyperbolic models or regions. However, profiles of the expansion 

I I scalar are less restrictive, with profile inversions necessarily taking place in 

elliptic models. We also examine radial profiles in special LTB configurations; 

CO closed elliptic models, models with a simultaneous big bang singularity, as well 

^ as a locally collapsing elliptic region surrounded by an expanding hyperbolic 

^^ background. The general analytic statements that we obtain allow for setting up 

^— H the right initial conditions to construct fully regular LTB models with any specific 

r^^ qualitative requirements for the profiles of all scalars and their time evolution. 

f"^ The results presented can be very useful in guiding future numerical work on 

• these models and in revising previous analytic work on all their applications. 

in 
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1. Introduction. 

m The well known spherically symmetric LTB dust models [1] are among the most 
useful known exact solutions of Einstein's equations. There is an extensive literature 
(see [21 13] for comprehensive reviews) using these models, as they allow for probing 
non-linear effects in inhomogeneous sources by means of analytic and/or tractable 
numeric solutions. In most applications LTB dust solutions lead to simple but effective 
toy models of cosmological inhomogeneities, as for example in the series of papers 
[H [3 El 13 IS|- The models are also useful in other theoretical contexts, such as 
censorship of singularities [H [TU] and even quantum gravity [TT]. More recently, 
LTB models have been extensively used in the widespread effort to explain cosmic 
dynamics without resorting to an elusive source like dark energy, either by fitting 
observations [l21[l3l[T4l[T5l[T6], or in the context of scalar averaging of inhomogeneities 
[TTI [T8l [T9l l20l [2Tj (see [221 ES] for comprehensive reviews on these applications) . For 
a novel theoretical approach to the dynamics of these models see [23] ■ 

Since large scale cosmic structure is dominated by large voids, the study and 
understanding of the spatial profiles of density and velocity in inhomogeneities is 
a relevant theoretical and practical issue that could provide important clues on the 
formation of these voids as part of the dynamics of cosmic sources. While LTB models 
provide an idealized description of cosmic inhomogeneities, a proper understanding 
(even qualitative) of radial profiles of scalars in these models still has a significant 
potential in astrophysical and cosmological implications in this context [51 [3] (see 
|25j for an update). In particular, knowing how these radial profiles evolve is 
extremely useful when applying to these models a formalism of scalar averaging of 
inhomogeneities, such as Buchert's formalism [501111]. 

An important theoretical issue in looking at density radial profiles is the possibility 
that a concavity change could occur under regular conditions (particularly around a 
symmetry center). This implies that a given initial density profile with a "clump" 
(over-density) form could evolve into a "void" (under-density) profile, or vice versa, 
a phenomenon that is usually denoted by a "clump to void" or "void to clump" 
density profile inversion. In a well known article Mustapha and Hellaby [IS] claimed 
to have furnished proof for the existence of this profile inversion, and to have found the 
conditions in which it occurs. However, since the conditions provided by these authors 
are too complicated (involving second and third order radial derivatives), their proof 
(their section 5) is based on restrictive assumptions, and thus is not applicable to 
generic LTB models. In the end, they discuss the observational implications of these 
profile inversion and furnish basically an empiric proof of their existence by looking 
at specific models and numeric examples where it occurs. 

Following a completely different methodology from that of Mustapha and Hellaby, 
we correct, extend and fully generalize their work in this article. By considering fully 
generic and regular LTB models, we examine the evolution of the radial profiles of 
the basic covariant scalars: rest mass density, p, spatial curvature, '^TZ, and Hubble 
expansion scalar, 8, aiming specifically at dealing in full generality with the following 
issues: (i) to provide a precise characterization of radial profiles as "clumps" or 
"voids", and (ii) to study in full rigor the conditions that allow for the existence 

f In order to motivate the reading of this long article, we have concentrated most of the background 
material in the Appendices and have written a proper summary of our results and their implications 
in the final section (section 11). Readers eager to know these results without going through the 
technical detail are advised to go directly to this section. 
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of profile inversions (assuming full regularity). As the reader can find out by going 
through the article (or by looking at the summary of results in section il), we do find 
fully analytic and rigorous conditions for the existence of clump/void radial profiles 
and their possible inversions for regular LTB models in general. 

An interesting and similar approach to the study of the evolution of radial profiles 
is found in the articles by Krasinski and Hellaby [31 [7] (see the review in [3]). These 
authors find the conditions for the existence of a unique LTB model, such that a given 
density or "velocity" profile at a given time constant hypersurface can be "mapped" to 
any other profile at a second hypersurface. While there is an obvious intersection and 
complementarity between the results of this article and those of [HIT], our methodology 
is different and (in our opinion) less restrictive, as we only require generic initial data 
to be specified at a single time constant hypersurface that can be considered a Cauchy 
surface in the context of an initial value problem. 

A summary of articles that have examined LTB models with void profiles in 
a perturbative and non-perturbative context can be found in the extensive review 
of this literature (previous to i994) given in section 3.1 of [2], and in the update in 
section 18.5 of [3]. Some authors [27l[28] have examined non-perturbative void profiles 
numerically in the context of structure formation scenarios (though, from our results 
in section 6, the elliptic models in [27] cannot be free from shell crossings). In some 
suggested configurations the void is modeled as an under-dense FLRW region, which 
is connected to a cosmic "background" through a section of a closed elliptic model 
[291 130]. In other articles [311 [32l [33] , the void is an LTB region containing a center, 
matched at a fixed comoving radius to another LTB region (the "envelope"), which in 
turn is matched to a FLRW "exterior". Although these models are rather artificial, 
their results agree with those of this article. 

The contents of the article is summarized in the remaining of this section. Section 
2 contains the bare basic background material needed to understand the article. The 
known analytic solutions of the field equations (given in Appendix A), which are 
determined by a preferred set of free functions (M, E, tbb)j have become the standard 
formulation in practically all theoretical, empiric and in numeric work on these models. 
While this standard parametrization is adequate and works in practice, we utilize an 
alternative set of covariant quasi-local scalars that are more suitable for our purpose 
(and that can be very convenient for numerical work, as well as a better understanding 
of several theoretical and practical issues) [M] 133 [M1I3Z1IM1I3S] ■ We also express these 
scalars and their fluctuations as quantities that are scaled with respect to their initial 
values in the framework of an initial value problem |38j . 

Extra background material, which is important as well but can be distracting if 
it appears in the main body of the article, has been placed in Appendices B, C and D. 
For example, the parametrization of the known analytic solutions in terms of the new 
variables and its initial value formulation (Appendix B), the Hellaby-Lake conditions 
[SHlliniHI] that guarantee an evolution free from shell crossings (Appendix C), and a 
discussion (Appendix D) of the relation between the radial coordinate and the proper 
radial length, which is the affine parameter along radial rays [3^. This is important 
since we are looking at the behavior of covariant scalars in the radial direction. 

We provide in section 3 a precise and rigorous characterization of the "clump" 
or "void" nature of a radial profile of a scalar A in terms of the sign of the gradient 
A' (monotonicity of A), which leads to the definition of a "turning value" (to be 
denoted "TV of A") as the value r = rtv where this gradient vanishes. In section 4 we 
examine the compatibility between initial clump/void profiles and the Hellaby-Lake 
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conditions, while the concept of "profile inversion" is defined in section 5 in terms of 
the vanishing of the fiuctuations defined in section 2, something which can occur with 
or without the existence of TV's. In section 6 we examine in detail density radial 
profiles and provide necessary and sufficient conditions for density profile inversions 
for parabolic, hyperbolic and open elliptic models or regions. We remark that the 
locus where the density fluctuation vanishes when there are density TV's corresponds 
to the "density wave" (see section 3.1 of [2]). The profiles and profile inversions for 
spatial curvature and the expansion scalar are respectively examined in sections 7 and 
8. The case of closed elliptic models is considered in section 9, as for these models the 
spherical topology of the hypersurfaces of constant time introduces a TV of the area 
distance R. In section 10 we look at radial profiles and profile inversions in special 
LTB configurations: simultaneous big-bang singularity, LTB models or regions not 
containing symmetry centers and configurations constructed by glueing LTB regions 
with "mixed" kinematics (elliptic, parabolic and hyperbolic regions, see [Uj), and 
specially the mixed configuration made by an elliptic "interior" region surrounded by 
a hyperbolic "exterior". Section 11 provides a summary and final discussion. 

Appendix E describes an alternative parametrization of the models in terms of 
quantities that can be related to the Omega and a Hubble parameters of a FLRW 
cosmology in the homogeneous limit. It is straightforward to translate the results of 
this article in terms of these parameters, which have been used in various articles 
[T^ [T71 [5S] describing cosmological inhomogeneities in various contexts (see also 
[38]). Appendix F provides the explicit expressions for the basic covariant scalars 
that characterize LTB models. 

2. LTB models, kinematic classes and a fluid flow time slicing. 

[U LTB dust models in their conventional variables are characterized by the following 
metric and field equations G"'' = Kpu°-u^ 

JDl2 

ds^ = -c^dt^ + dr"^ + R^{de'^ + sin^ Odip^), (1) 

R' -^-f+E, (2) 

2M' = KpR^R', (3) 

where u'^ — 6q, p is the rest-mass density, n = SttG/c^, E = Eir), M = M{r), R = 
R{ct,r), while R' = dR/dr and R = u^V^i? = dR/d{ct). 

The analytic solutions of Q are given in Appendix A as parametric forms 
involving R, M, E and a third function tbbir) (the "big bang time"). It is common 
usage in the literature (see [2 El ES HI]) to classify these solutions in "kinematic 
equivalence classes" given by the sign of E, which determines the existence of a zero 
of R^. Since E = E{r), the sign of this function can be, either the same in the 
full range of r, in which case we have LTB models of a given kinematic class, or it 
can change sign in specific ranges of r, defining LTB models with regions of different 

§ This section provides the minimal background material to make this article as self-contained as 
possible. The reader is advised to consult references |34ll37ll38l for details on the quasi-local scalar 
representation of LTB models. Necessary background material is also found in the Appendices. 
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kinematic class (see [IT]). These kinematic classes are 

E = 0, Parabolic models or regions (4a) 

E > 0, Hyperbolic models or regions (46) 

E < 0, Elliptic models or regions (4c) 



where the equal sign in (46) and (4c I holds only in a symmetry center 



1 — r'" T-j-\i r^ mn/j™ ' ^ ' 



2.1. Covariant quasi-local scalars. 

The normal geodesic 4-velocity in (IT]) defines a natural time slicing in which the space 
slices are the 3-dimensional Riemannian hypersurfaces ^T[t], orthogonal to u°, with 
metric hab = UaUb + gab, and marked by arbitrary constant values of t. All radial rays 
in a given '^T[t] are diffeomorphic to each other and to the real line M (or to continuous 
subsets of R if there are symmetry centers) . Hence, we will consider every LTB scalar 
function as equivalent, under the time slicing given by u°, to a one parameter family 
of real valued functions A[t] : K -J> K so that A[t](r) = A{t, r). 

An LTB model is uniquely determined by a representation of local fluid flow 
covariant scalars A (see Appendix F), though alternative scalar representations can 
be defined. For every scalar function A in LTB models admitting (at least) a symmetry 
center (r = 0), we define its quasi-local dual Aq as the family of real valued functions 
Aq[t\ :M+ ^M given by [ii] 

j^AFdVp _ gAR^R'dx 
/J'j'dVp ^ j^ K^R'dx 

where the integration is along arbirtary slices ^T[i],with J^ = (1 + -E)^/^, dVp — 
^/habdrdOd-Lp, and we are using the notation L ...dec = J Zq ...dx. 

The definition of quasi-local scalars ([5| leads in a natural manner to an initial 
value parametrization of LTB models, so that all quantities can be scaled in terms 
of their value at a fiducial (or "initial") slice '^Ti = "^Titi], where t = ti is arbitrary. 
Hence, the subindex i will denote henceforth "initial value functions" , which will be 
understood to be scalar functions evaluated at t — t.^ . This scaling of quantities to '^Ti 
suggests introducing 

L^^, (6) 

SO that dependence on R becomes dependence on L as a sort of dimensionless scale 
factor. By introducing the following definitions (which simplify notation) 

2m=-p, 2mq = -pq, ^="5"' « = ~6~' 3"' "^ ^ ~3~' ^' 

it is straightforward to obtain from ([2]), ([s]), ( |F.1[ ), ([s]), ([6| and ([t]) the quasi-local 
density, spatial curvature and expansion scalar 
ruq, M 

h F 

' L2 i?2' ^^> 

nl = JJ^ = —^=2mq-kq= ^^^3 "' , (10) 

II While we have assumed that the integral in 1151 is evaluated for a fully regular ^T[t], it is 
straightforward to generalize this definition to hypersurfaces that intersect singularities and to models 
without symmetry centers (see Appendix A-5 of [38]). In order to simplify the notation we will omit 
henceforth the symbol [t] attached to scalar functions unless it is needed for clarity. 
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so that the conventional free ftinctions M and E are expressible in terms of initial 
value functions by 



M = mq,Rl, 



E ■ 






(11) 



2.2. Fluctuations. 



The quasi-local scalars defined by ([5]) comply with the following properties 



A', = {A,y 

A{r)-Mr) 



3R' 



[A~A,], 



1 



A' R^ dx. 



Hence, given the pair of scalars {A, Aq], we define their relative fluctuations as 



5(A) 



A-Aq 
A„ 



A'/A 



q/^'-Q 



1 



3R'/R Aq{r)R3{r 



A' R^ dx, 



(12a) 
(126) 

(13) 



where we used (12a) and (12&). The local density, spatial curv ature (Ricci scalar of 
the '^T[i]) and expansion scalar 8 — Vau"', given by ([3| and (F.l ), are then expressible 
in terms of the fluctuations (13 1 or their gradients (from (12a|): 

m' 

(14a) 



m — rrif, 



rv — fv/-! 



ri — rLn 



1 + ,5(") 

l + JW 

l + <5(«) 



« ' 3R'/R' 

K 

^ "^ 3R'/R- 



Ho 



n' 



iR'/W 



(146) 



(14c) 



The remaining local fluid flow scalars in (F.2) and (F.3|, associated with the shear 
and electric Weyl tensors, follow as 



By introducing the scale factor 

r.4^^1 



F — -m (5^™) — ^— 



L'/L 



(15) 



(16) 



R'i/Ri R'i/Ri 

the following scaling laws for the local density and spatial curvature follow readily 
from (|3]), ([9|, ^ and (|f1]) as: 

[l + ^(")]^i^, (17a) 



iii-qi 



L3r' 



fc = ^[i + ^^^-^] = ^ 



r-i 



3(1 



icC^)^ 



(176) 



Hence, by comparing (17a| and (176) with (14a| and (146), we obtain the following 
scaling laws for the fluctuations (5^™^ and S''^' 

1 + (5(") = ^ , (18a) 



2 I j(.0 _ 2/3 + Jf ) 

3 r ' 



{iim 



(186) 
(18c) 
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while the scahng law for (J*^^-* follows from (13), (14c) and (15) 



2^(„) ^ 2m, (5(") - kg (je^) _ 2m,, ,5(") - kg.LS^'^^ 



jIiTTLq Kq 



ZfTlqi kqi h 



d"^) 



i-r]-fc,,L[5f) + |(i-r)] 



(19) 



[2mqi - kqiL] r 

which allows us to obtain scaling laws for the local expansion scalar, "H, and the scalar 
functions S and £. 



2.3. LTB metric and curvature singularities. 

It is important to remark that the scalars m,, T-Lq, kq and their fluctuations are 
covariant objects, as Af, E, R, R = u'^VaR are invariants in spherically spacetimes 
[45] , Given (^, ^ and ([l6]), the LTB metric ^ takes the form 

r^ Rf dr^ 



ds^ = -c\\t^ 



L' 



R^ (d6|2 



sm 



(20) 



The Friedman-like equation ^ now takes the form (10), whose solutions are 
summarized in Appendix B and are equivalent to those of ([2|) in Appendix A. However, 
the solutions of (10) are expressed in terms of L and the initial value functions 



rriqi, kqi, Ri, with the "big bang time" t^^(r) given in terms of the latter functions. The 
analytic solutions of ( 10 ) in Appendix B allow us to compute the following analytic 



form for F by implicit radial derivation of (B.l ), (B.5) and (B.IO), and then using (13) 
and L'/L = (1 — T)R'jRi to eliminate rn' j^,k'^ and L' in terms of 5™ 
The result is 



(5?^ and F 



• Parabolic models or regions. 



F = 1 + 5i™^ - ^' 



(m) 



L3/2- 



• Hyperbolic and elliptic models or regions. 



?(™) 



T=i+3{sr' -sr'){i- 



(fe)^ 



Hq 



qt 



-'dUqCit^t,) 



^irn) _ 3 (fe) 



(21) 



(22) 



where Tig and Tiqi follow from (10), while c{t — ti) is given by (B.5) and (jB.lO) 



By inserting (21) or (22) (depending on the kinematic class) into (|17a)-(17&). 



(18a)-(186) and (19) leads to closed analytic expressions for all these scaling laws 



in terms of L and initial value functions in the context of the initial value approach 



described in previous subsection. The functional forms of F and Hq in (10) and (21) 



([22| allow us to obtain analytic expressions for all the scaling laws of covariant scalars 
(as functions of L and initial value functions) . 

The scaling laws (l8])-( 19 ) clearly indicate the existence of two possible curvature 



singularities whose coordinate locus is 

L{t, r) — 0, central singularity 

F(i, r) = shell crossing singularity. 



(23a) 
(236) 



so that for reasonable initial value functions (bounded and continuous), all scalars 

m, k, % can also 



A„ 



m. 



q, kg, T-Lg diverge as L — )■ 0, whereas local scalars A 



diverge if F — >• (even if L > 0). 
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Notice that if F > 0, then aU scalars A and Aq only diverge at the central 
singularity L = 0, which is an intrinsice feature of LTB models. However, there is 
a shell crossing singularity if F — )■ for L > 0, so that all the relative fluctuations 
Jf"^' and local scalars A diverge while their quasi-local duals Aq remain bounded (with 
Aq ^0). This is an obviously unphysical effect of shell crossings that must be avoided. 
We will denote by "regular LTB models" all configurations for which shell crossing 
singularities are absent, thus complying with 

F > V {ct, r) such that L>0 (24) 

In order to test this regularity condition we need to check the sign in the forms of F in 



(21 1 and (22). However, as proven by Hellaby and Lake [301 SI] i these conditions 
can be specified in terms of M, E, R and their gradients, which in terms of our 
parameters translates into being specified in terms of initial value functions rriqi, kqi 
and their initial fluctuations 5^ , 5i . The Hellaby-Lake conditions are summarized 
in Appendix C separately for parabolic, hyperbolic and elliptic models or regions. 

3. Qualitative features of radial profiles of scalars. 

3.1. The radial coordinate. 

The study of radial profiles of covariant scalars (local or quasi-local or fluctuations) 
is basically the study of their behavior along radial rays in the hypersurfaces '^T[i] 
orthogonal to u". Since the radial coordinate has no intrinsic covariant meaning, 
while the radial rays are spacelike geodesies of the LTB metric p^, we should use 



their affine parameter (proper radial length) £ given by (D.l) to probe the radial 
behavior of covariant scalars. [J] Though, unfortunately, the dependence of scalars on 
i is extremely difficult to evaluate, even qualitatively. However, as long as we assume 



that ( 24 ) holds (shell crossings are absent) and r is a well defined coordinate so that Ri 



and kqi comply with (D.2 ), the proper length ^ is a monotonously increasing function of 
r. As a consequence, dependence on i becomes qualitatively analogous to dependence 
on r (see |39|). and thus, all results given in terms of r regarding features of radial 
profiles of scalars (monotonicity, TV's or asymptotics) will be qualitatively analogous 
and fully equivalent to results given in terms of i. We will henceforth assume that 



24|) and ( |D.2[ ) hold. 

The fact that the radial coordinate in (fTl) or (20) admits an arbitrary rescaling 



r = r{r) defines a 'radial coordinate gauge freedom' that allows us to use any function 
of r as a radial coordinate (for example setting M as radial coordinate as in [H |21 [El [Ij 



and many other articles). In particular, the metric form (20) suggests using this 
coordinate gauge freedom to prescribe a functional form for the initial value function 
Ri{r) [3lJ [38l [39] . However, the form of Ri is not completely arbitrary: it is restricted 
by the topology of the space slices '^T[t] and by regularity conditions (see Appendix 
D and also [33] and Appendix A3 of [55]1. 

3.2. Monotonicity and characterization of radial profiles as "clumps" or "voids". 

Besides their role as covariant non-linear perturbations |37l I38j . the relative 
fluctuations 6^ ' convey important non-local information on the radial profile of 

^ See Appendix D for a discussion on the regularity of the proper radial length and the radial 
coordinate. 
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scalars A. Consider the radial domain (i constant) 

^^{r) = {x\0<x<r}, (25) 

where the lower bound r = marks a symmetry center. We shall assume henceforth 
that R' > holds for every domain of this form. We have then: 

Lemma 1 . Let A be a smooth integrable scalar function defined in i?c('')- If 
A is monotonous => Aq is monotonous, and the following results hold for its 



quasi-local dual Aq and relative fluctuation (5^"^' defined by ([5J) and ( 13 1 

• sign(A - Aq) = sign(Aq(5(^)) = sign(A'). 

• li A^O, then sign((5(^)) = sign(A7A). 

Corollary. The following sign relations hold: 

• li A>0, then Aq > A when A' < and Aq < A when A' > 0. 

• li A<0, then Aq < A when A' < and Aq > A when A' > 0. 



Proof. These results follow directly from (12a), (12&) and (13). The converse 
statements are not true in general (see the comment to Lemma 3). 
Comment . We will consider (unless stated otherwise) domains of the form 
^c(^) whose lower bound is a symmetry center. However, Lemma 1 and 
its corollary are valid for any regular radial domain, in particular domains 
not containing a symmetry center. Such domains arise in either one of the 
following situations: 

• When looking at a given region in compound configurations made by 
elliptic, hyperbolic and parabolic regions [H] (see section 10). 

• LTB models that do not admit symmetry centers (see section 10). 

• Radial domains of hypersurfaces '^T[t] that intersect a central singularity 
marked by ct — c<bb(^) or ct — ct^^n{r) in the {ct, r) plane (see Appendix 
A-5 of [35]). 

In all these cases the domain must be a suitable modification or restriction of 
'&c{i'). In the case of the intersection of ^T[t] with a singularity, the domain 
can given by i9(r) = {x \ r^^^^ < x < r}, where r^i^^ is the value of the radial 
coordinate marking the intersection of ct^,^^(r) or cf<,oi,(r) with a given ^T[t]. 
Evidently, d{r) does not contain a symmetry center. 

Comment: clump and void profiles. Lemma 1 indicates how, for a monotonous radial 



profile of A in a domain (251, the sign of (J'"*' determines if this profile has the form 
of a "clump" or "void". ['^I Consider a non-negative quantity such as 2m = Kp/3, a 
density clump (m' < 0) in a domain i9c(r) around a symmetry center of an arbitrary 
hypersurface '^T[t] will be characterized by (5^™-' < 0, whereas for a void (m' > 0) we 
have (5'^'"^ > 0. If we consider the case of negative spatial curvature (fc < 0), then the 
signs of fc' invert: fc' > for a clump and fc' < for a void, but in either case the 
sign of 5^^^ remains the same as that for positive quantities: (J^™-* < for a clump 
and (5*^™^ > for a void. In case there is a zero of Aq, making S^'^'' diverge, then the 
"clump" or "void" nature of a radial profile should be characterized by the sign of 
A ~ Aq — Aq5'^^\ Since the results of this lemma are valid for an arbitrary t, they 
can be very useful for characterizing the features of initial conditions at a given t ^ ti 
and their evolution for t ^ ti. 

+ This characterization becomes more nuanced when dealing with domains not containing a center, 
and also when R' = holds without violating regularity conditions (in closed elliptic models in which 
the ^Tlt] have S'^ topology). See sections 9 and 10. 
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3.3. Turning values (TV's) 

If the radial profile of a scalar in a given '^T[t] is not monotonous, then there exists 
(at least) one turning value of A (a "TV of A"), which is a value rtv G ^c{r) such 
that ^'(ftv) = 0. We discuss below the conditions for the existence of a TV. Let A 



be a smooth integrable scalar function defined in the domain 'dc{i') in (25) along an 
arbitrary ^T[i], then 

Lemma 2 . The necessary and sufficient condition for the existence of a TV 
of Aq in a hypersurface '^T[t] is 

<5(^)(rf,rtv) = 0, (26) 



Proof. The result follows directly from ( 13 1 given the assumption that R' > 
holds in every dc{r). Since Aq = Aq{t,r), then the value of rtv will (in 
general) change for different ^T[t] (see figure Q. 

Corollary 1. If r^v & iSdr) marks a TV of Aq, then Aq{rty) ~ A (rtv). 



Corollary 2. From section 4 of [39]: if (D.2) holds, then 3 r = rtv such that 



^'('"tv) = ^ 3 i = itv such that dAjai = at ^ = £tv, where i is proper 
length (|D.1[) along the ^T[t]. 



Comment 1. When A is the density, the locus of (261 in the {ct,r) plane is 
the so-called "density wave" mentioned in section 3.1 of [2j. 

Comment 2. The case in which R' — occurs under regular conditions ('^T[i] 
homeomorphic to §^) is discussed in section 13. Notice that 5^"^^ -> ±oo with 
r > marks a zero of Aq, not (in general) a TV of Aq (see Appendix A4 of 

EH]). 

Lemma 3 . The existence of a TV of Aq in 'dc{r) is a sufRcient condition for 
the existence of a TV of A. 



Proof. Considering the corollary of the previous lemma and using (126) we 
have at rtv 

v4(rtv) - A,(rtv) = /''" A'{x)^^dx = 0. (27) 

Since R'^{x) > for x > and _R'^(rtv) are positive, then for this integral to 
vanish A' must have a zero in the range < a; < rtv. 

Comment. The converse of this lemma is not true: a TV of A does not imply 
a TV of Aq in every domain. This is connected to the fact that a monotonous 
Aq does not imply (in general) a monotonous A (the converse of Lemma 1). 
Consider a domain ■!?c(''tv)- The TV of A occurs at a value x = ri < rtv, 
therefore for values ri < a; < rtv the gradient A' has already changed signs 
whereas A' has not. However, if we consider domains with r > rtv and there 
is a TV of A at X = ri and i?' > in i}c{r), then for sufficiently large r there 
will always be a TV of Aq at x = rtv > fi . 

It is clear, as a consequence of the lemmas proven above, that knowledge on the 
qualitative form of the radial profile of Aq in a given '^T[t], i.e. monotonicity, TV's 
and radial asymptotics (see next section), immediately translates into knowledge on 
the radial profile of its dual scalar A. As shown in section 2, the quasi-local scalars 



Evolution of radial profiles in regular LTB dust models. 



11 



rUq, kq, T-Lq satisfy scaling laws that are less complicated than those of their local 
counterparts m, k, H. Hence, it is more practical to examine first the radial profiles 
of the scalars Aq in order to use the lemmas of this section to infer the profiles of the 
A. Finally, we remark that monotonicity or existence of TV's can change from one 
^T[t] to the other. 

4. Initial profiles compatible with absence of shell crossings. 



It is worthwhile looking carefully at the types of profiles (clumps or voids) co mp atible 
with initial value functions rriqi, kqi (or T-Lqi, Vii in Appendix E) that satisfy (24), and 



so assure a time evolution free from shell crossing singularities. For this purpose, 



we need to examine the allowed signs of 5\ and 5\ involved in the Hellaby-Lake 



(m) 



e(fe) 



conditions (C.6), (C.7) and (C.8| 



It is evident (from (131 and Lemma 1) that clumps or void initial profiles directly 



correlate with the sigs of the gradients tti'j, fc' j, but it is not so clear how this type 
of profiles can be prescribed by means of the conventional free functions M and E. 
To appreciate the relation between (5^™-', (5*^'^-' and their initial values 5, 



(m) .(fe) 



to 



the functions M and E and their gradients, we use the scaling laws (18a I and (186) 
together with (13), (Isl) and (|9|, leading to the general form of the relation (C.3) 



M' 


3R 


M 


R^ 


E' 


3R[ 


E 


R, 



r[i- 

2 
3 






Ri 
Ri 



[1 + ^. 

2 
3 



Ml 



?(fc) 



(28a) 
(286) 



From these relations, we can convey the characterization of clump {5l < 0) or void 
{Si > 0) initial profiles to the following constraints on M and E and their gradients: 



Initial 



Initial k 



M'/M 
iR'^/Ri 
E'/E 
2RJR, 



< 1 (clump) 



< 1 (clump) 



M'/M 
3R'jRi 

E'/E 



> 1 (void). 



> 1 (void). 



(29a) 



(296) 



where E > and E < hold, respectively, for hyperbolic and elliptic models or 
regions. 



The parab olic c ase is trivial. Since fulfillment of (|24|) only requires the Hellaby- 

must have 



s-(m) 



Lake condition (C.6): —1 < S^ < 0, the initial value functions ' 



^qi or rtqi 



clump profiles, so that M must comply with the restriction in the clump case in (29o) 
The hyperbolic and elliptic cases require further examination. 



4-1- Hyperbolic models or regions 



The Hellaby-Lake conditions (C.7) do not rule out a clump or void profile of rriqi and 



kqi. The first two conditions simply place a lower negative bound for S^ and 5, 



(m) 



e(fc) 



hence to get further information we look at i^,^ < from (C.5) rewritten as 



-/3 5f^ <0, 



where 

1 j~LqiC\Z^ 



a 






[2- 



11/2 



arccosh(l 



3/2 



(30) 



(31a) 
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P = 7;^q^c{t^ - 4b) - 1 = X + 



3 



3[2 



11/2 



arccosh(l 



2x: 



3/2 



(316) 
(31c) 



f^qi 



Since a > and /3 > for every choice of initial value functions (either rriqi, \kqi\ 



'Hqi, Vti), then (30) necessarily excludes the possibility that 51 > 0, 5l > 0, and so 



(m) 



c(fe) 



initial void profiles for density and curvature are incompatible with the fulfillment of 
( |24|. O n the other hand, sufficient (but not necessary F] conditions for the fulfillment 
of (24) follow by choosing initial density and curvature profiles of clumps: 



1 < ^i™^ < 0, 

s-(m) 



< a(") < 0. 



(32) 



Still, initial density voids 5l > are compatible with (24) provided initial curvature 



is a clump: 5\ < 0. The opposite situation: density clump and curvature void is also 
compatible with (24). In fact, by looking at the functional forms in (31a) and (31&), 

0.327), with: 



we have a ~ (3 for xt « 4.11 (or fi. 

a> P for < a;^ < 4.11 or < Cl, < 0.327, 
a< P for Xi > 4.11 or 0.327 < (li < 1, 



(33 a) 
(336) 



where Cli = 2mqi/'H'l^ ~ 2/(2 + Xi) (see Appendix E). Hence, other sufficient (but 



not necessary) conditions for ( 24 ) involving initial voids of m and k can be given by 



selecting x-i to be restricted to specific ranges: 



Am) 
i 
Am) 



> 0, (5f ^ < 

< 0, 5^'"^ > 



with 
with 



Am) 
i 






(k)i 



(m). 



for 
for 



< f^, < 0.327, 

0.327 < a < 1, 



(34 a) 
(346) 



In general, given a choice of rriqi, \kqi\, a very practical way to test the fulfillment of 



(24) follows by rewriting the Hellaby-Lake conditions (C.7) as 



> 



vria 



Ri 



tZa 



> 



21^ 
Ri 



to', x,a{xi) + fc', ^{x,) < 0, 



"-gi S 0, 



(35) 



E'>0, as in (CD 



where a and /3 are given by (31a) and (316), and Ri can always be specified as a 
coordinate choice (for example , as i n (D.3)). The constraint (35) can be given in 
terms of fli, Hqi by means of (E.18). In terms of M and E, we have M' > and 



It can be argued that choosing the free functions M, E and t^^, directly from the 



Hellaby-Lake conditions, as in (C.I), is easier than testing (30) or (35). However, the 



compensation for going through these more complicated constraints lies in the fact 
that, for some specific problems, it can be more practical and/or intuitive to select 
initial Cli and Hqi, related to the Omega and Hubble factors (see Appendix E), or an 
initial density and curvature profile, than functional forms for M, E and ^t- 



4-2. Elliptic models and regions. 



It is evident that the Hellaby-Lake conditions (C.8) are far more restrictive than in 



the hyperbolic case (C.7). In order to examine the type of initial profiles (clumps or 



This corrects an error in |46i . where these conditions were regarded as necessary and sufficient 
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voids) compatible with (24 1, we rewrite the first two conditions in (C.8) as 
with «!, /3i and a2, P2 given by 



ai 



/3i = 



a2 



2 (2 — .Ti)^/^ arccos(l — Xi) 
^ 3/2 ' 

6 — a;i 3(2 — Xi)-'^/^ arccos(l — Xi) 



2xi 



3/2 



(2- 



U/2 



2a;: 



[27r — arccos(l — Xi)] 



3(2- 



3/2 
)i/2r27r- 



arccos(l — cc^)] 



ZXi 



2x, 



3/2 



(36 a) 
(366) 

(37a) 
(376) 

(38 a) 
(386) 



where < Xi < 2 (which is equivalent to i7,; > 1). Since ai, /3i and a^, P2 are all 
positive for every choice of initial value functions tti, 
these conditions into a single inequality: 

/32 



qi, kqi or r^i, Hqi, wc Can combine 



a2 ' "1 



< ^<5?' 



where ^>^>1>^>0. (39) 

2 a2 ai 



Since this inequality is valid for all initial value functions rriqi^ kqi or Hqi^ Cli, then it 
clearly implies the following necessary (but not sufficient) condition for an evolution 
free from shell crossings complying with ( 24 ) 



s;(™) 



<o, 



<0. 



(40) 

Hence, as opposed to the hyperbolic case, the Hellaby-Lake conditions ( C.8[ ) are 
incompatible with initial voids of density and curvature. The functions M and E must 
then comply with the clump restrictions in (29a)-(296). Notice that (C.8 1 imposes 
the restriction 5^ > — 1, but places no lower bounds on a negative (5^ 



:{k) 



The constraints (36a)-(366) and (39) are too complicated for guessi ng th e form 

q^ or -Kqi, Cl, to fulfill (las I. 



It is 



of the appropriate initial value functions rriqi, k, 

more practical (see [IS]) to test any given ansatz for these functions first on the 



necessary conditions (C.lll, (C.IO) and (40), and once verifying that these conditions 



are satisfied, the ansatzes should be tested on (39), which can be written in a form 



similar to ( 35 ) 






< 



< 






(41) 
qi a2 rriqi kqi ai 

where ai, a2, Pi, P2 are given by (37a)-(376) and ( |38a )-(386) and we used (13) to 
express (5^ , 5l in terms of m! i < and k'^ < (though niqi and kqi are necessarily 
non-negative). This condition (which is equivalent to the Hellaby-Lake condit ion 
(C.2)) can be rephrased in terms of H'i, Ct'^ by means of the relation (E.18) in 
Appendix E. 
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Figure 1. Profile inversions with turning values. Panel (a) displays 
the {ct,r) plane associated with a radial domain | |25| . The region in which 
^C™) (^ct^ r) < is shaded, so that its upper boundary marks the surface 
5("»)(ct, r) = (upper thick curve). The locus of the initial singularity ctbb(^) is 
the lower thick curve. The black dots A, B,C mark the points where (5'™' = 
intersects r = r2, r^, r^ where there is a turning value of niq for ct = ci2> cia, ct^. 
Panel (b) displays qualitative plots of the radial profiles of niq normalized to the 
central values mq^ = mq{ct,0) for the same values of ct and r shown in panel 

(a). Notice that an initial "clump" profile was selected, so that &^ < and 
^„i "^ hold in all the radial domain, but for times ct > etc this initial profile 
has evolved into a void profile with (S'™) > and m' > 0. Notice (from Lemma 
3) that turning values of ruq imply turning values of its local dual scalar m. The 
situation described in this figure for rrtq is applicable to turning values and profile 
inversions of other scalars like kq or Hq in hyperbolic models (as S^'^' diverges as 
"Hq — > in elliptic models, see section 8). Notice that the surface 5'^"^>(t,r) = 
is the so— called "density wave" [2]. 



5. Profile inversions and turning values. 

In reference [Hj Mustapha and Hellaby addressed the problem of whether a density 
profile of the type of a clump(void) may be "inverted" into that of a void(clump) as 
a regular LTB model evolves in time. These authors examined the concavity of the 
radial density profile around the center of symmetry, but their results involve quite 
cumbersome second and third order radial derivatives of R and the free functions 
M, E, ct^^i,. Their proof of the existence of this inversion (their section 5) is based on 
restricting assumptions and is not fully general. In the remaining of this article we 
extend and complement these results and show by using Lemmas 2 and 3 that the 
issue of a "profile inversion" can be addressed rigorously and in a much more intuitive 
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cf 


1. (a) 


c4 














8'"' > 




6r=o 






Ctt 




8'"' < 


cf^ 






rt^ 






Cti 




''•■bb 





mjmqc 




Figure 2. Profile inversions witliout turning values. Panel (a) displays 
the {ct,r) plane associated with a radial domain ( |25| . The region in which 
S('^>(^ct,r) < is shaded, so that its upper boundary marks the simultaneous 
surface S^"^'{r) = 0, which corresponds to ct = cti (upper thick line segment). 
The locus of the initial singularity ctbh{r) is the lower thick curve. Panel (b) 
displays qualitative plots of the radial profiles of rtiq normalized to the central 
values niqc = mq{ct, 0) for the same values of ct and r shown in panel (a). Notice 



that 5, 



(m) 



implies m' ■ = so that rUq 



marking ct = cti the locus of 



the profile inversion, as mq has clump profiles for ct < cti and void profiles for 
ct > cti, without density TV's in any value of t. Notice (from Lemma 3) that 
turning values of niq imply turning values of its local dual scalar rn. The same 
picture is applicable to profile inversions of kq or Hq. 



way, not only for density but also for spatial curvature and the expansion. 

In order to address the "clump into void" question in precise and rigorous terms. 



we consider a quasi-local scalar function Aq defined in a radial range ( 25 1 in a regular 



LTB model complying with condition (24), so that R' > {'^T[t] homemorphic to 



see section 9 for the case of ^T[t] homeomorphic to S"^). 



Definition: Profile inversion . The scalar Aq undergoes a profile inversion in 
a domain ^c(^) if there exist two values t — ti, ^2 with ti < t2, such that for 
all X e -dcir) 



sign 



Si^) 



for t < ti ^ sign 



SiA) 



for t > t2, 



(42) 



The "clump into void" inversion of Aq corresponds to the particular case 
J(^' < for t<ti and (5^-^) > for t > ta, (43) 
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where the equal signs above hold only at the center. The opposite inversion 
(void into clump) is defined by the opposite signs. 

Comment . It is useful to identify ti with ti, so that the sign of the initial value function 
6l conveys the initial clump or void nature of the profile that becomes inverted for 
some ^2 > ti. However, ti < ti < t2 is also possible, which means selecting Aqi whose 
profile is not monotonous. 

Lemma 4 . The existence of a profile inversion of Aq implies a profile inver- 
sion of A. 

Proof . Following Lemma 1 and its corollary, a clump (void) of Aq in a given 
domain ^c(?') implies a clump (void) of A. Hence, the result follows. It is 
important to remark that for practical purposes it may be far easier to ex- 
amine the profile inversion on Aq, so that the existence of a profile inversion 
of A follows from Lemmas 1 and 3. 

Lemma 5 . A sufficient (but not necessary) condition for a profile inversion of 
Aq in a domain t?c(?') is the existence of TV's marked by a; = rtv[t] in every 
^T[t] with t € {ti < t < ^2} and with either one of rtv[ti] = or rtv[i2] = 
holding. 

Proof . Consider a clump into void inversion for ^^ > and the case 
''tv[i2] = (the void into clump inversion and the cases rtv[ii] = or neg- 
ative Aq are analogous). Since there is no TV before ii, then A' < and 
5^^) < hold for all x G -ddr) for t < ti. For each tj e {ti < t < ^2} there 
is a TV at a; = rtv[ij] so that 5^^\ctj,rtv[j]) — 0, and so we have: A'^ < 
ior < X < rtv[tj] and ^^ > for rtv[tj] < x < r, so that rtv[tj] ^ rtv[tfc] 
for every tj ^ tk. As a consequence, ([26]) is not marked by a surface of 
simultaneity in the (ci, r) plane. Since there are no more TV's for i > <2 in 
'(?(r), then A' > and (5^^^ > hold and the clump profile has been inverted. 
This process is illustrated by figure [T] 

Corollary 1. Lemma 5 also provides a sufficient condition for a profile inver- 
sion of A. This follows from Lemmas 2 and 3, as the sufficient condition for 



the emergence of a turning value {A' = 0) at some r — rtv is simply ( 26 1 



Corollary 2. If the initial value function Aqi is selected so that A' ^ = S^ ' = 
for r = rtv[ti] € "ddr) and there exist t > ti such that (5^'*'' 7^ for all t > i and 
r > 0, then there is a profile inversion of Aq in the domain -ddrtviti]) C '&c{t)- 
From Lemmas 3 and 4, there will also be a profile inversion of A. This is a 
restricted form of Lemma 5 in which there is no information on the sign of 
J(^) for t<t,. 

Comment . The choice of an initial value function Aqi with a non-monotonous 
profile (there exists a TV such that 5l =0 holds for some r — rtv[ij] G 
dc{r)) is not sufficient for having a profile inversion, since 5^ ' may keep 
the same sign for all i > t^ in the domain "d c{^tv[ii\) ■ In general, further 
information on 5^ > is needed given a non-monotonous profile of Aqi. 
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We have provided sufficient conditions for a profile inversion of A based on the 
existence of a: 
without TV's. 



existence of an inversion of Aq and TV's, however, such inversions can also occur 



Lemma 6 . If there are no TV's oi Ag, a sufficient (not necessary) condition 
for a profile inversion at a fixed t = tg follows if ( 26 1 is marked by a surface 
of simultaneity t — t^ 

(5(^)(cio,a;)=0 V x e ?9c(r), 



(44) 



so that (5^"^^ has opposite signs for t < tQ and for t > Iq. 



Proof . Consider the case of a clump into void inversion and Aq > 0. Since 
there are no TV's, then d*-"^^ < in i!}c{r) for all t < to and (5^^^ > in 'dc{r) 
for all t > to- From Lemmas 1-3, the sign of 5^"^-' determines the sign of A' 
hence (since there are no TV's) a profile inversion occurs at t = to- This 
situation is illustrated by figure [2J Notice that S^'^') = at t = to does not 
imply that 6'^^'^ vanishes at t ^ to. 

Given the definition of a a profile inversion and its relation with the existence of TV's 
in the radial profile of scalars, we need to examine now the conditions for S^^^ — for 
the main quasi-local scalars Aq — mq, kq, Jiq. As we show in the following sections 
these conditions are easier to analyze and more intuitive that the concavity test around 
the center considered in [26]. 



6. Radial density profiles and their inversions. 

Following Lemma 1, the qualitative clump or void form of the radial profile of 
niq = {K/3)pq is determined by the sign of S^""^^ . From Lemmas 1 and 2 this information 
readily translates into general statements on the profile of local dens ity m = (k/3)/9. 
To examine the sign of 5'-™-' we obtain this function directly from (18a I, which by 



using ( 22 ) and ( 22 ) leads to 



j(m) _ "t ^ 



(m) 



rL3/2 ^^^g 



(m)^ 



g{m} ^ :^ 



-(m) 



2 ' 



^3/2 „ j|™) ' 
HqCit^t,) 



parabolic 



:(»") 



Ak) 



T~Lq 
Ha 



hyperbolic and elliptic 



(45 a) 



(45 &) 



where F, "H^ and Hqi are given by (10), (21), (22) and (E.5), and we assume that 
initial value functions 5^ \ 5\ and l-Lqi have been selected so that the Hellaby-Lake 



conditions ( |a7| and ( [081 h old (F^O) 



Substituting F from (21) and ( |22[ ) into (45a) and (456), and bearing in mind that 
->■ oo as i — )■ ibb (or L — )• 0, see Appendix B), we find that in this limit 



^ -1, (46) 

holds for parabolic, hyperbolic and elliptic models and for whatever initial value 
functions that we might choose at t = t^. The same limiting value holds as i — >■ t^^a in 
the collapsing singularity in elliptic models (where "Hq — >■ — oo as L — > 0, see Appendix 
B). From (18a), the limiting behavior (46) necessarily implies that F — )■ oo as L — >■ 0, 
as in general S™ > — 1. We also assume domains ddr) given by (25) in which kqi is 
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either zero or does not change sign in the full domain, which means that 'i?c('') defines 
a "pure" parabolic, hyperbolic or elliptic region containing a center (possibly part of a 
"pure" model) in which 6^ is bounded (mixed elliptic/hyperbolic configurations are 
examined in section 10). 



6.1. Profile inversion: "Density clumps into voids" revisited. 



A density profile inversion of mq in a regular LTB model complying with (24) is 



characterized by ( 43 1 for A 



From Lemma 4, ( 43 ) also implies a profile inversion 



for the local density to. Considering domains of the form i9c(^), sufficient conditions 
for the profile inversion of mq can arise if 

• There exist density turning values (to,' = and S^"^^ = 0) at r = rtv[t] for slices 
^T[t] marked hy ti < t < t2 with ti < ti (Lemma 5, figure [l]). 

• If there are no TV's, then there must exist t = to such that S^"^\cto, x) = and 
the sign of (5^™^ for t > to is the opposite to that for t < to (Lemma 6, figure llj). 



We examine the conditions for the fulfillment of these two possibilities. 

Consider the case of regular parabolic models. It is evident from ( |C.6 ) and (45a) 
that (5^™^ < must hold in any radial domain ■&cir) in these models (with the equal 
sign holding only at the center). Hence a density profile inversion is not possible: 
their density profiles are monotonously decreasing for all t so that initial parabolic 
clumps remain clumps for all the time evolution. For domains in hyperbolic or elliptic 



configurations the Hellaby-Lake conditions (C.7) and (C.8) do not exclude a priori 



the possible emergence of a density profile inversion, but checking if conditions ( 43 ) 
or 



( [44]) can be fu lfilled requires further examination. In practice, we need to examine 
how (5^™^ in (456) evolves in time given the existence (or lack) of turning values. This 



process is illustrated by figures [T] and [2] 



6.2. Hyperbolic models or regions. 

6.2.1. Inversions with TV's We assume henceforth that initial value functions have 



been selected so that (24) and thus, the Hellaby-Lake conditions (C.7), hold for all 
t. For arbitrary fixed (and finite) r and large times t/t^ ^ 1 we have from (10): 
L « j/cqil^/^, which leads to L sa \kqi\'^^^it/ti) ^ 1. Hence the large time asymptotic 
regime follows from asymptotic series for L ^ 1: 



T~Lq ^ 
Hqi (2 



X 



1/2 



Hqc(t-t,)«l 



2 hiL 
Xi L 



+ 0{L-'), 



0{L-^), 



(47a) 



(476) 



where c(t — ti) follows from (B.5). These asymptotic forms imply that Hq/Hqi and 
T-Lqc{t — ti) in (22) and (45 6 1 respectively decrease and increase monotonically between 
and 1. The asymptotic for of (5'™^ for t/ti ^ 1 is 

jf"^ - (3/2)5f ) r 



l + (3/2)jf) 

where we used ([22|, ( |456[ ) and ( |47a[ )-(|476| 
result: 



2 Ini 
Xi L 



0(L- 



(48) 



Therefore, we have the following general 
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Lemma 7 . If (5^ " =^ 0, 6l 7^ comply with the Hellaby-Lake conditions 



(C.7|, a necessary and sufficient condition for the existence of a "clump to 



void" density profile inversion in a regular hyperbolic model is 

Si^) _ 3^(^-) > 0, (49) 

where the equal sign holds only at the center. 
Proof . In the asymptotic regime t/ti ^ 1 the form of the leading term of 



^(™) in (48) becomes independent of Xi (or m^j, kqi) for sufficiently large L, 
hence (49) implies (J*^™-* > for all r in this limit and the converse statement 
is also true. On the other hand, 5'^™-' — > — 1 as L — > for all x G i^dr) is 
necessary and sufficient if the Hellaby-Lake conditions to hold. Hence, ^'■"'■' 



must change sign at least once if (49) holds and the converse statement is 
true (see figure fl] for a graphic depiction) . 

Corollary. A "void to clump" inversion without shell crossings is not possible. 

Proof . A "void into clump" scenario requires the combination S" > and 
(5('n) < 0, which would only have a chance to occur if 5\ > 0, but this initial 



condition together with (5^ > violates (30), and so is incompatible with 



r > (and consequently, with the Hellaby-Lake conditions). 

Comment: no profile inversion. The converse statement of Lemma 7 implies 



that the reverse condition, 5^ — (3/2)(5,- ' < (plus (C.7)), indicates an 
shell-crossing free evolution without a density profile inversion: the initial 
density clump profile is preserved (see Table 1). Notice that a void profile for 
the whole time evolution is not possible because near the initial singularity (as 
t -^ tbb) we have necessarily S^"^") — > —1 for every choice of initial conditions. 
However, the pure void evolution is possible for models with a simultaneous 
big bang (t^^ = 0) for which (5^™) -^ as i — > ttt (see section 10). 

We examine now the consequences of Lemma 7 for different assumptions on (5^ and 
Si . If we choose an initial density with a clump profile (— 1 < i5j < so that 
5j = — |<5j I) and desire an evolution with a clump to void evolution, then the 
profile inversion must occur for times t > ti, irrespectively of the existence of density 



TV's. In this case, condition (49) implies 



Jf)<0, |5f)|>?|eY (50) 

If we assume an initial density void (5^ > 0, together with L > 0, then a profile 
inversion must have already happened for previous times tbb < i < ti because 
J(™) ^ _i as t -^ ibb (with ot without density TV's). Also, if 5^™^ > and F > 
hold, then (48) implies an asymptotic void profile ^'■'"^ > if (5^ < 0. The allowed 



qualitative evolution of the radial profiles of density and spatial curvature are d epic ted 



in Table 1 for all possible sign combinations of (5- ™' and S- ' compatible with (|24 ) . 



6.2.2. Inversions without TV's While Lemma 7 applies to either one of the possible 
situations described by Lemmas 5 and 6, the clump into void inversion without density 
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TV's is compatible with less restrictive forms of the initial value functions. Since t — to 
is arbitrary, is is convenient to examine this case by assuming that t^ — ti, leading to 

Lemma 8 . If there are no density TV's along the hypersurfaces ^T[t] of a 
regular hyperbolic model, a sufficient condition for a clump to void density 
profile inversion in 'dci'f) is 



Proof, lis. 



(m; 



<5i"^ = 0, 
= 0, then (5(^ 



i<* 



(fe) 



<0. 



j(™) — 



u 



(k) r 



1 



Hn 






in (456) takes the form 

. (3/2) <5f) 



l + {3/2)S. 



(k)-' 



(51) 



(52) 



where we substituted S^ — in the expression for T in ( 22 ) and used only 



the leading term in the asymptotic expansions (47a)-(476). Evidently, S^ 



in (|52|) implies the conditions on S^ specified in (51). Since S^" = and 



S("^)> as t — > oo, and there are no density TV's, then S^™'' > holds 
for all t > ti. On the other hand, since (5^™^ — >■ — 1 as i — >■ itb (or L — > 0), 
then the lack of density TV's implies that (5^™-' < for all t < ti. Hence, a 
clump to void (and not a void to clump) density profile inversion must occur 
at t = ti (see figure [2]). 



6.2.3. Initial conditions that yield inversions. It is important to remark that initial 
conditions complying with a clump into void density profile inversion are not strange 
nor unusual, and are perfectly compatible with the Hellaby-Lake regularity conditions 



(C.7). Condition (491 can be expressed in terms of gradients of initial value functions 



by the equivalent inequalities: 



3 kqi 
9 h 
3E' 



M 2 E 



iii-qi 
M' 



> 0, with: m^ < 0, fc^^ > 0, kg, < 0, 



> 0, with: 



„ M' 3i?' 



°<5^f 



m 



while condition (51) is equivalent to 
= 0, fc;,; > 0, k, 

Ri 



qi 
M' 
Id 



<0, 



°<f^f 



(53 a) 
(536) 

(54a) 
(546) 



Notice that these conditions are compatible with the constraints for an initial clump 



in regular models given in (29al and (296|. Also, we can rewrite (53a|-(536) and 



(54a|-(546) in terms of the radial coordinate by means of the coordinate gauge (D.3) 



6.3. Elliptic models or regions. 

The following general result holds for elliptic models or regions in which R' > holds 
everywhere (no TV of R) : 

Lemma 9 . There are no density profile inversions (with or without density 
TV's) in any domain "ddf) in a regular elliptic model or region lacking a TV 
of R. The proof of the is given below. 



Evolution of radial profiles in regular LTB dust models. 



21 



Corollary. Density has a clump profile for all the evolution time and for every 
domain 'dc{r), hence —1 < (5'™^ < (clump profile) must hold everywhere. 



If Lemma 9 holds, the proof follows directly from (18a| and (C.8). 



Proof of Lemma 9 . We examine first the conditions for a profile inversion when there 



are density TV's. Regularity implies that F > holds for all times, so that (55a) 



(55c I hold. Hence, in order to examine the conditions for a density profile inversion 



in these configurations, it is useful to rewrite conditions (C.lOl, (39) and (40), which 
are necessary and necessary and sufficient for F > 0, as 



Sim) 



(m) 






- Pi ' 



% r^^-l^rl, 


(55a) 


>0 ^ |5f^|>3l<5f)|. 


(556) 


l"''\<7^'\Sn with: "^^^>"^>L 
/3ia2 " Pia2 /3i 


(55c) 



where ai, a2. Pi, 
rewrite (|45&|) as 



'2 are given by (37a)-(376| and (38a)-(386). Considering (55a), we 



where 



-l'^. 



(m). 



3C 



1-3C 



C=|5f''|$-|,5, 



(fc) 



*<3' 



T~Lo 



$ = l--/-H,c(i-i,), 

itqi 



n-qi ^ 



u), 



(56) 

(57a) 
(576) 



From the expressions above we can readily examine the behavior of f^'™' at maximal 
expansion, t = t,„^^ (or L = L,-„„,, = 2/xi) where Hq = 0, as well as near the collapse 
t — >■ icon where Hq 



—00. The result is 

('^i-(2/3)|4"'^| 



Sim) ^ _^\_S 

1 



3\5., 



(m). 



^(" 



-1, 



t 



toon. 



t^t„ 



(58a) 
(586) 



Comparing (58a) with the necessary condition (556 1, it is evident that F > implies 
that (5(™) < at i = i„„,. Since (5^™) < at i = t„,, and as t -^ t,„„, the only 
possibility for S^"^^ > to occur is for some intermediate times tt < t < t„,^ and/or 
tma,x < t < icoii- We now assume that (5'^™-' > and F > both hold for these times 



and verify if this assumption is consistent with (55a )-(|55c[ ) 
If <5(" 



> and F > both hold we have from ( 56 



1-5, 



Wi 



< 



3$ 



1|5W| 



<\s. 



(fe)i 



1-|'5> 



(m). 



(59) 



^■q, ' ^ ' ' ^ ' 3* 

an inequality that can be examined qualitatively by means of the analytic expressions 
(|B.10|) and (|E.12|). In the range t, <X< ^".- the term (3$ - l)/(3*) grows 



monotonically from —00 to 2/3, hence (59) implies that \5l \ < {2/3)\Sl \ holds 



irC^)! 



)C(™)| 



in all this range, which is in contradiction with the necessary condition (556 1 and thus 
implies that 5'^™^ > cannot hold. 
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Now consider the range t^^^ < t < i^„„, in which the term (3$ — l)/(3^) grows 
monotonicaUy from 2/3 to a maximal value as L — > (or t — >■ i^„i,) given by 



3*(0,a;,) ~ /32{x^) 



<1, 



>-qt 



2{n, - 1) 



(60) 



where < Xi < 2 and a2, 1^2 are given by (38a|-(38&), while Cli is defined in 
Appendix E. Since now (59) is not in contradiction with (|556 | because it allows for 
\5\ I > 2/3|(5j I, we need to see if it is compatible with (55c). Combining the latter 



and ( 59 ) leads to the inequality 

a2(a;i)|„(„i) 



\5r\<w'\< 



(fc)i 



3$(L,x, 



\5.. 



(m). 



(61) 



h{x^)' » ' ' » ' 3*(i,a;,) 

whic h sh ould hold for all Xi and L. Si nce (pO| is valid for all i nitial value functions, 
then (61) is clearly inconsistent, and so (59) is inconsistent with (55cl and so (5'^™-' > 



cannot hold in this range of t. 

For a profile inversion without TV's we need to impose the condition 5"^ = 0, 
as in (51) (see Lemma 8). From (B.12) and (C.5) and (C.9| two of the Hellaby-Lake 



conditions in (C.8) now take the form 



ct' 



ZR[/Ri 
ZR[/Ri 



(fc) 



3Ze(a;,)(2-x,)i/2 



2x. 



3/2 



'^f^c(i. 



tbb) + 



ctL 



3R[/Ri 



<0, 



>0. 



(62a) 



(626) 



Since the term in square brackets in (62a I is negative for all the allowed range 
< Xi < 2, fulfillment of this regularity condition requires S^ > 0, but then condition 



(62 6 1 will not be satisfied. Therefore, there are no regular elliptic models complying 



with the conditions for a profile inversion without TV's. 



7. Radial profiles of spatial curvature. 

Just as with the density, the qualitative behavior of radial profiles of kg — ^TZq/Q and 
k = ^TZ/fS are determined by the sign of 5^^\ which takes the following form for both 
hyperbolic and elliptic models (with '^T[t\ ^ M'^): 



j(fe) 



2 
f 



I (^5(™) _ 3^(.)^ [^^,(, _ ,^) _ 1] + (5(™) _ s^)^ 






(63) 



where F, Uq and c(t-ti) aie given by pij| ), ( [21) ), ([22]), pl5[ ) and ( [B.10| , and fulfillment 
of regularity conditions (C.7) and (C.8) is assumed. Notice that 6^'^^ = holds for all 



domains 'dc{r) in parabolic models and regions. 

As in the previous section, we now substitute F in ( 22 ) into ( 63 ) . Bearing in mind 
that V-a 



oo as i — >■ ih 



j(fe) ^ _ 



(or L 
2 
3' 



0), we find that in this limit 



(64) 



holds for hyperbolic and elliptic models and for whatever initial value functions that 
we might choose at t = i^. The same limiting value holds as i — >■ t^^^ in the collapsing 
singularity in elliptic models (where T-Lq 
behavior (64) necessarily implies that F 



— oo as i — > 0). From (186), the limiting 
oo as L — > 0. 
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Case 


^t-) 


,(.) 


r(m) 3c(fe) 


Sira) 


S(k) 


Comments 


i 


<o 


<o 


>o 


>o 


<o 


No profile inversion of k. 

Profile inversion 

of TO at t > ii 


ii 


<o 


<o 


<o 


<o 


>o 


No profile inversion of to. 

Profile inversion 

of fc at t > ij 


iii 


<o 


>o 


< (nee) 


<o 


>o 


No profile inversion of to. 

Profile inversion 

of fc at t < ii 


iv 


>o 


<o 


> (nee) 


>o 


< 


No profile inversion of k. 

Profile inversion 

of TO at i < ii 



Table 1. Initial conditions and radial profiles of density and spatial 
curvature in hyperbolic models. The table presents a summary of the 
qualitative evolution of the radial profiles of m and k for all possib le co mbinations 



of 5. j^ and 5 ^ 0, compatible with regularity conditions (C.7l. The fifth 



and sixth columns display the asymptotic sign of <5('"' and (5(*). The term "nee" 
means that 5^ — (3/2)S^ ' necessarily has the indicated sign for these specific 
sign choices of <5 and 5 . 



The profile inversion of spatial curvature was not considered by Mustapha and 
Hellaby, but this issue can also be examined along following the same arguments of the 
previous section. Table 1 summarizes the possible qualitative evolution of k and kg for 
all combinations of initial value functions compatible with the regularity conditions 



(C.7) and (C. 



Hyperbolic models or regions. 



Bearing in mind the asymptotic expansions (47a)~(476) and T given in (22 1, the 
form of (JC^) in (|63| for i > 1 is 



1 + (3/2)jf ^ 



2 InL 

Xi Li 



+ 0{L-') 



(65) 



Hence, 6^'^^ — > as L -^ oo (or i — )■ cx)), but its sign is determined by the oppo- 
site condition to the sign of (J*^™-' , which means that the features of the radial profile 
evolution, including conditions for a profile inversion of k, are given by the complemen- 
tary situations to those of to.. As a consequence, the sign combinations of (5j , S^ 
and (5j — (3/2)(5j- that yield a given sign for d^"^^ , will yield the opposite sign for 5^'^^ . 

Elliptic models or regions. 



As in the previous section, we rewrite conditions (C.IO), (39) and (40) as in 



(55a)-(55c), transforming (k33| into the following form similar to (56) 

{k 



-1^1 



2C 



1-3C 



(66) 



Evolution of radial profiles in regular LTB dust models. 24 

where C is given by (56) and (57a)-( |576|. Si nce $'■'''> -^ -2/3 as i -> ttb and t -> i,„u, 
whereas S^'^^ tends to the same hmit (58a| as i — > t„„^, then foUowing the same 



steps as in case of (5'"*'' in the previous section we can prove that S^''^ < holds 
for all the evolution range tbb < i < ^con- Therefore, just as with density profiles, 
spatial curvature profiles in regular elliptic models whose hypersurfaces ^T[t] have M.^ 
topology have the form of clumps for all t. 



8. Profiles of the expansion scalar 



From (10), (14c), the scaling law (19) and (E.6)-(E.8) we can readily infer the sign 
of (5'-^^ from what we have found in previous sections about the signs of ^'■'"^ and 
j(fc), Por parabolic models we have 26^'^^ — (5^"^ hence regularity conditions require 
gCH) < Q (^j g^ g^ clump radial profile for H — 8/3 and Hq — 6g/3) for all the evolution 
time t > ibb- 

For both hyperbolic and elliptic models or regions we have 2mq/'Hg = 1^ — >■ 1 as 
L ^ {i.e. as t -^ i^b and t — >■ ic„„ in elliptic models). Hence, in this limit we have 
§{-H) ^ ^(m)/2, and so 

<5(«)^-i. (67) 

We examine now the sign of S^'^' for domains 'i9c('") in hyperbolic and elliptic config- 
urations. 



Hyperbohc models or regions. 



From ( |E.7| and ([19|, we have < fi < 1 and 51 -> for i > 1 (or t -> (xd). 
Hence J^^-* — ?► (5'''''/2 in this limit, and so (5^^) has the same asymptotic sign as (5''^' 



(see Table 1). Bearing in mind (67), there must be a profile inversion of H and Hq 



<c("i) 



for the cases (ii) and (iii) in Table 1 for which S^/"' - (3/2)5f'^ < and 5^''^ > holds 
asymptotically: 



Case ii 
Case iii 



girn) < Q Ak) < Q 



\sl'^\<l\sl"'^\ 



(68 a) 

(686) 
From (19), we have S^"'' < in case (ii) where both Si"''' and 5\'^' are non-positive, 



r(m) 



< 0, 5f ^ > 0, 

e(W) 



(m) 



M 



hence in this case the profile inversion must occur for some t > ti. However, for case 
(iii) it is not evident if this inversion happens before or after t = ti. We examine this 
case below. 



AH) 



5?^^ < 



The sign conditions of S^ for the case (iii) in (68&) can be given by 
2 



c(fc) 



-i < 



?(«) 



>0, 



1-5. 



C™)! ~ X, 



profile inversion at t > ti, 



profile inversion at t < ti, 



(69 a) 



(696) 



r(™)| 



e(m) 



where \S\ | < 1 (because S^ < 0). We need to verify the compatibility between 



(69a)-(696) and the regularity condition (30), which for the sign combination (686) 

a{xi) 



takes the form 



Ak) 



\5 



(m). 



< 



{^^)' 



(70) 
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where a{xi) and (3{xi) are given by ( |31a[ ) and (31&). By looking at these expressions, 
we can readily see that 2/xi < a{xi)//3{xi) holds for the full range of Xi, which can 
be expresse d as < Cli < 1 where fli = 2/{2 + Xi) (see Appendix E). Hence, all initial 
conditions ( 69a ) associated with S^ < will be compatible with (30 1, whereas initial 
conditions (|696|) associated with S^ > need to satisfy the extra constraint 

a{xi) 



< 



(k) 



(m) 



< 



P{X^)' 



(71) 



As shown by these expressions, the sign of 61 (and thus the time of the profile 
inversion of Ji) is not obvious for the case (iii), and so it must be found by testing the 



fulfillment of (69a)-(696) and (71) for any given choice of initial value functions. 



A regular evolution without profile inversion of H is also possible. The necessary 
(not sufficient) condition for this is given by cases (i) and (iv) in Table 1. However, 
it is not possible to provide a condition for this evolution without going to particular 
cases or specific restrictions of the initial value functions. 



Elliptic models or regions. 



Near the initial and collapsing singularities (ibb and tcou) we have S^^'' -^ ^1/2, 
but since Hq — >■ and f2 — >■ oo as i — )■ t^^x, then 6^^^ diverges in this limit. We 
examine then the radial profiles of the expansion by looking at the sign of 



''«*"" -«-««^5Fi«- 



2mq (5(") - 

mo 



kg J^*^) 



(72) 



where we used (13), (E.1)-(E.2) and (E.5). Notice, from (15l, that this quantity is 



the n egati ve of the scalar S associated with the shear tensor (see F.3)). From (19), 
(22 1, (456) and (63), it is straightforward to show that the leading term oi T-LqS'-^'lor 



L„ 



IS 



y a(«) -- 




LL. 


14 


vith- '''- Yjf") ^ 




^'^*^- 3RyR. - [^^ 2 


whereas as L -> (as i — > ttt an 


At 



ct' 



,(m) 



,(.)) 



SR'jRi 
ct' 



nj^^^ 



' L3/2 



> t^oii) the leading terms are 

TOO, 



(73) 



(74) 



(75) 



where c{t 



(± 



tbb) is give n in terms of initial value functions by (B.12|-(B.13) or 



of Hq in the limits L 



Tiq takes the opposite sign 



( |E.16| . The T sign in (JTSJ) indicates that "H,^'^' = V. 

0, so that < "H < Hg as t ^- tbb (when H, Hg — )• cx)) and 
teoii (when v., T-iq -^ — cx)). This is consistent with (5^^' -^ ^1/2 



■H, < -H < as i 
in both limits. 

Evidently, irrespectively of the choice of initial value functions and of the sign 
of ct[^^^/R!^, the fact that T-LqS^"^^ evolves from — oo to oo impHes that there must 
necessarily be a change of sign of T-L' and so a profile inversion of T-L is an inherent 
feature of open elliptic models that expand and collapse. We have at least the following 
two possibilities: 
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If tLx/^» > 0> then -Hq(5(^) > as t ^ i„ 
changed sign at some t < t^^^. 

If iLxZ-^i ^ 0; tlien ■HqfjC") < has not changed sign as t 
the change of sign of V.qS'-^'' must occur at some t > t„^^^, 
as i — >■ t^„,,. 



and so this quantity has necessarily 



— >■ i„^^. In this case, 
so that y.qS^'^^ -)■ oo 



In order to test these possibiUties, we need to examine sign of ct'^^^/ R[, which is 
easier to do by looking at the function ct^^^{r) given by (B.13) and by considering the 



possible radial behavior of Xi = kqi/ruqi and j/. 



:,l/2. 



^ai. ^' 



in open elliptic models 
Considering (|D.2|) and using the coordinate gauge ( |D.3[ ), kqi is restricted by 

„2 



1 



E = I- kqir- 



Ro>0, 



<Xi <2, 
which implies that (i) the limit k, 



qi 



(76a) 
(766) 

-^ fco = const. > as r — !> oo is not possible, so 
that kqi must decay radially like r^^ or faster, and (ii) if niqi — >■ then ruqi must 
decay as fast kqi, though rUqi — >• toq = const. > is possible. This yields the following 
possibilities for the asymptotic limit of Xi (see [39] ) 

mo, Xi -> 0, j/j -> 0, (77a) 

(776) 



kq, -^ 0, 
kqi -^ 0, 



^0, 



Xi~^ xq< 2, yi -^ 0. 



Notice that (77a) implies Xi < 2, as x — ?> 2 is only possible if assuming (11 h). If the 



asymptotic behavior is given by (77a) or (77&) with xq < 2, then it — Zf,{xi) remains 



positive and bounded, and thus ct 

TT 

> OO. 

= kn/mo 



ct„ 



Cti 



If we consider (776) with xq 
can assume that both m 
around e = we obtain 



(r) in (B.13) takes the radial asymptotic form 

(78) 

2, then x, w 2 



e were e ^ 
can assume that both mqi and kqi decay as r""*" with 2 < 7 < 3 (see 



ct„ 



cti 



Zi{xi)] 



.7/2 



2^/2^ 



ct. 



/Too 



0{r '>'), as we 
I). Expanding 



(79) 



As a consequence of (78), if rriqi and fcg^ behave asymptotically as in (77a) or (776) 
with xo < 2, then we must necessarily have ct'_^^^ > at least in the asymptotic radial 
range, though ct'^^^ < might still occur for values of r closer to the center. However, 



in the case (776) with xq = 2 it is impossible to assert the sign of ci^^^ without making 
further assumptions on rUqi and kqi. Radial profiles and profile inversions of Hq are 
displayed by figure [3] for the case ct[^^^ > 0. The profiles for the case ct'^^^ = would 
be analogous to those of figure [2] depicting a clump to void profile inversion without 
TV's, as in this case Hq = is simultaneous. 



9. "Closed" elliptic models. 



So far we have only considered radial domains i9c{r) defined by (25 ), containing a single 
symmetry center and complying with R' > 0. The upper bound of these domains can 
be extended to encompass full "open" elliptic models whose hypersurfaces ^T{t) are 
homeomorphic to M.^. If the hypersurfaces ^T{t) are homeomorphic to S'^ {i.e. a 
"closed" model), then there is a second symmetry center r = r^, so that the full radial 
range is {0 < r < rd, and there is necessarily a TV of i? at ?■ = rtv < re such that 
R'irtv) = occurs in regular conditions. 
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Figure 3. Profile inversions of the expansion scalar in elliptic models. 

Panel (a) displays the {ct,r) plane associated with a radial domain ( |25| . The 
region in which Hq {ct, r) > is shaded, so that its upper boundary marlcs the 
maximal expansion surface t = tmaxir), which coincides with 'Hq{ct,r) = and is 
not (in general) simultaneous. The locii of the initial big bang and final collapse 
singularities, ctbb(^) and ctcoii(''), are the lower and upper thick curves. Panel 
(b) displays qualitative plots of the radial profiles of Hg for the same values of 
ct and r shown in panel (a). Notice that the change of sign of Hn and the fact 



that "Ho 



itoo at the singularities strongly constrain these radial profiles. Since 



"^^max ^ 0, then (see section 8) the profile inversion has already taken place for 
ct4 < cimax (where the equality only holds at r = 0). From Lemma 3, the profiles 
of the local expansion H qualitatively behave as those of Hq, but the surface 
H = does not coincide with ct = ct„iax(^). In closed elliptic models Hq and H 
have a common TV with _R at r = rtv (not shown here). For these models the 
curves in both panels are qualitatively the same in the range < r < rtv, and are 
mirror images reflected by a vertical axis for rtv < r < re, where r^ marks the 
second symmetry center (see section 9). 



All domains 'ddf) with r < rtv in closed elliptic models are wholly equivalent 
to elliptic domains examined so far, hence the results proven in previous lemmas 
concerning elliptic models and regions remain valid. In particular, from Lemma 9 
and the results of sections 6 and 7, the density and spatial curvature profiles in these 
domains must be monotonous and of a clump form for all t. However, domains i^df) 
with r > rtv will contain a TV of R, and thus the results of the previous lemmas must 
be re-considered. We examine in the remaining of this section the properties of radial 
profiles in domains with r > rtv 



From the regularity condition (D.2|, a well defined proper radial length (D.ll 
requires the ratio R' /J- with T — \/l + E to be continuous (at least C°) at r = rtv. 
Since i?'(rtv) = -^(rtv) = 0, with R"{rtv) < 0, J"'(rtv) < 0, then from I'Hoppital rule 
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we have at the hmit r — ^ rtv 

R' i?"(rtv) 



>0, 



(80) 



J^ -F'(rtv) 

which imphes that rtv is also a TV of the quasi-local volume element: V'dr = 
AnR'^R'dr, but not of the proper volume element: V^dr = 47r(_R^i?'/J^)dr. The 
spherical topology of the ^T[t] places a strong constraint on the radial profiles of 
scalars. We prove now the following 

Lemma 10 . Let A be any one of the covariant local scalars {m = np/Z, k = 
^Tl/Q, H = 6/3} associated with domain ^dr) with r > rtv in a closed 
elliptic model. Then 



i?'(rtv) - 



^'(rtv) = 0, 



(81) 



The converse statement of (81) is false: as we proved in previous lemmas, 



TV's of TO, fc, H may exist under the assumption that R' > holds. 



Proof . By considering (12 a I and (80), it is evident that for all Aq 
i?'(rtv) = => 4(rtv) = 0. 



(82) 



To prove (81 ) we assume that rtv is a TV of R and examine the behavior of 
A and Aq around r = rtv + e for |e| ^ 1- Considering (from (82)) that a TV 
of R implies ^^(rtv) = 0, we have at leading orders: 

Aq{rt^ + e) « Ag(rtv) + -A','(rtv)e^ 

i?(rtv + e) « i?(rtv) + ^R" {rt.)e^ , 

(83) 



Since, in general, A(rtv) 7^ 0, by applying (12a) to A and Aq and at leading 
orders we get 



^(r-tv + e) - A{rt^) = Aq{rt^ + e) - A<,(rtv) + 



A'qR 



3R' 



i"tv+e 



A'qR 

3i?' 



(84) 



Using ( 83 ) to evaluate this expression at leading orders and dividing both 



sides by e we get 

A(rtv + e) - A(rtv) _ A, (rtv + e) - Ag(rtv) , 1 ^„ 

6 



+ 7?<(^tv)e. 



(85) 



e e 

Taking in both sides the limit as e — >■ we get A' (rtv) = ^n(^tv) — 0, which 
is the desired result. 



Corollary 1. A TV of i? at r = rtv implies a TV of M, E, ttb, icon and t,„^^ 
at r = rtv, but not a TV of J^. Bearing in mind that a TV of R in closed 
models holds for all t, it implies f?^(rtv) = 0, the proof then follows directly 
from ( [C^ , dCj] ), (p8^-([28&| and ^. On the other hand, a TV of £; does 
not imply a TV of J^, which satisfies 1 ~ T'^ = kqR^ . If we differentiate 
both sides of this expression and consider that the regularity condition (D.2 ) 



requires that i?'(rtv) — and J^(rtv) = hold, we find that R'{rtv) does not 
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imply J"'(rtv) = 0. 
Corollary 2. A TV of i? at r = rtv does not imply a change of sign of 5^^^ at 



r = rtv The proof follows directly from the definition of 5'^'^^ in (13 1 which 
involves the quotient of A' and R' . At r — rtv the sign of A' and A' in 
the range < r < rtv changes (Lemma 10) to its opposite sign in the range 
rtv < r < re, but the sign of i?' changes at rtv in the same manner, thus 
keeping the sign of 6^ ' in < r < rtv equal for the whole radial domain. 

As a direct consequence of Corollary 2, 5^"^'' < and (5''^'-' < must hold everywhere 
in all regular closed elliptic models. Since domains i)cif) with r < rtv are equivalent 
to domains considered by Lemma 9, then 5*^™^ < and (5*^'^-' < must hold in these 
domains, but since the TV of i? at r = rtv does not change the signs of (S*-™^ and J*^'^^ 
then J^™) < and 5'^'"'^ < must hold also for every domain with r > rtv 

The common TV of scalars at the TV of R (Lemma 10) marked by rtv is an 
inherent feature of the S"^ topology of the '^T[t], and as such it is independent of the 
time evolution. Once we select an initial ^T[ii] with this topology (which requires 
choosing Ri such that i?^(rtv) = 0), then rtv will be a common TV of the initial value 
functions rriqi, kqi and (from Lemma 10) of m^, ki. The TV at rtv is then common 
to m, niq, k, kq, Ji, T-Lq. While there is no profile inversion of the density and spatial 
curvature because 5^™-' < and S'^^^ < must hold for all t, there is a local minimum 
of m and rriq at rtv, and so we have a sort of "topological density void" in the form of 
an under-dense region marked by coordinate values around r = rtv, between the two 
symmetry centers. This topological void is a sort of thick annular shell of under-dense 
dust layers concentric to both of the symmetry centers, and so it is of fundamentally 
different nature from the "normal" density void examined in previous sections, which 
occupies a closed compact under-dense region around a symmetry center. However, 
the most important difference is the fact that the normal density void is not a feature 
occurring at a fixed r = rtv, but an evolved feature at changing values r — rtv[i] (the 
"density wave" [5]) that depends on specific generic initial conditions at i = i^, hence 
it always arises as the produce of a 'clump to void' profile inversion. 

Since regular closed elliptic models exist that contain an arbitrary number of TV's 



of R (as long as (D.2| holds at each TV), any of these configurations constitutes a 
complex pattern of topological clumps and voids in which the density profile changes 
from an over-density to an under-density and back at the TV's. However, as in the 
case of a single TV of R, these regions are concentric under-dense or over-dense thick 
spherical shells, not compact regions, and are fixed features (not evolving "density 
waves") that follow from the choice of an initial ^T[ii] in which initial value functions 
Ri, rriqi, kqi have an arbitrary number of common TV's. This inherent feature of 
closed elliptic models to produce TV's of density has been used in the literature to 
generate void models by means of compound configurations constructed by suitably 
matching regions of closed elliptic models and sections of FLRW dust models ^29, _30] . 
Since the voids and clumps in these configurations do not follow from generic initial 
conditions, they are wholly artificial and will not be considered any further. 

Regarding the profile of the expansion scalars H and Jiq, we take into 
consideration that the region < r < rtv of any closed model is qualitatively analogous 
to a region of an open model containing a symmetry center. Therefore, the results 
of section 10 hold for such regions, which then necessarily have a profile inversion 
associated with a TV of Hq, which depending on the sign of t'^^^, will be either in the 
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expanding or collapsing phase. From Lemma 3, there is necessary a profile inversion 
from a TV of "H as well. Since H and Hq must have a TV at r = r^v, then radial 
profiles in the region rtv "£ r < re are mirror images (reflected by r — rtv) of the 
profiles in the region < r < rtv (as those depicted by figure Isl) . 



10. Special configurations. 

10.1. Models with a simultaneous big bang. 

Regular hyperbolic and elliptic LTB models exist for which the initial singularity 
occurs in a single singular ^T[t], hence i^b = ^o = const. < t^. We prove in this 
section that density and spatial curvature profile inversions are not possible for these 
configurations. 



Considering (|B.6|), (|B.7|), (|B.ll|), (|B.12|), (|C.5|), a simultaneous big bang implies 
follow 

c{ti - to) 



the following constraints among the initial value functions 



c{ti - to) 



Zh{x^) 


gi^) 
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P 


Vi 
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ZeiXi) 


gim) 


- ^j^*^) 


0< 



>o, 



Vi 



«! 



"1 



<1, 



hyperbolic, 



elliptic, 



(86 a) 
(866) 



where Zh, Zf. are given by (B.4) and (B.9|, while a, /3, ai, /3i were defined by (30) 



and (36a). The constraints (86a)-(866) are clearly incompatible with (51), hence there 



are no density profile inversions witho ut a T V (t urnin g va lue). 

The functions F, ^('"' and (jC^' in (Eil, (EH&I and (63| take the form 



F = 1 



-m 



rtqi 



- "iUq C{t - to) 



^im) _ 3 ..; 



(87) 



.(m) 



r(m) 



2 ' 



"Hqi c{ti — to) 



[Hq^ C{t, - to) - 1] 



The Hellaby-Lake conditions (|C.7| and (|C.8) reduce to [35] 

hyperbolic, 



<5i"' > 



<5i™^ > 



-1, 



-1, 



5?^ > 



sr(™) 



^^f ^ > 0, 



elliptic. 



(89) 

(90 a) 
(906) 



We look at the hyperbolic and elliptic cases separately below. 

• Hyperbolic models. We have 2/3 < Hqc{t ~ to) < I with Hqc{t — to) — > 2/3 and 
Hqc{t — to) — > 1, respectively, as L — ?> and L — > oo. Hence, if the Hellaby-Lake 

conditions hold (F > 0) the sign of (J^™' is the same as sign of Sl"^' — {3/2)Sl , 
but from (86a) the sign of S^' must be opposite to the sign of S^ ' , and also we 
have 

Am) _ 3 (fe) ^ _ //3 3^^ ^(fe) 

2* U 2, 



(91) 



If ^1™^ < and (5f ^ > 0, then 5|™^ - (3/2)(5^ < follows from (jgij) and so 
Sim) < Q^ jf §{rn) ^ q ^^^ g(k} ^ q^ ^^^^ ^ ^.^^^^ ^(„,) ^ ^ ^^ obtain 

analogous results for spatial curvature, since the sign of S*-^^ is opposite to the 
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sign of S\" — (3/2)(5| < 0, so that 6^''^ has the same sign of S^ ' . Therefore, 
there are no density or spatial curvature profile inversions, though a density void 
profile is possible for all t > to, as (5^™^ — >■ as L — > at the simultaneous big 
bang t = to. 

• Elliptic models. Since Hqc{t — to) = at i,„ax, we have < Hqc{t — to) < 2/3 in 
the expanding phase and 'Hqc{t — to) < in the collapsing phase. The constraint 
(866 1 implies that 5) must have the same sign as 5^ , hence the Hellaby- 



Lake conditions (|906|) can only be fulfilled if —1 < 5-"^ < and (5- < hold. 



Following the same reasoning as in the hyperbolic case, it is easy to show that 
^("i) < and iJ^*^) < must then hold for all t. Therefore, there are no density 
and spatial curvature profile inversions, and (as opposed to the hyperbolic case) 
density void profiles are not possible in regular configurations. 

10.2. "Mixed" models in which kg changes sign and models lacking a center. 

Regular LTB models can be constructed by glueing regions of hyperbolic, parabolic 
and elliptic types, either by matching the parts at fixed comoving radii (as in 
[Sn [551 1551 m] ) , or simply by prescribing as part of the initial conditions a function 
kqi that changes sign along its radial domain. It is important to remark that matching 
conditions for this type of kinematically mixed regions allow for step discontinuities 
in the radial direction and on radial gradients of scalars. Thus, since S^"^\ ^'■'^■' 
depend on the radial gradients m'^ and fc^, these relative fluctuations could be only 
piecewise continuous in these constructions (unless the sections are glued with the 
extra requirement that radial gradients are also continuous). 

The results proved in the lemmas in this article remain valid when applied 
separately to each region of any compound configuration made by glueing any 
combination of parabolic, hyperbolic or elliptic regions complying with regularity 
conditions. This is obvious for any region (for whatever sign of kq) containing a 
center r = 0, in which ^'''^•'(t, 0) = must hold for all t. However, the lemmas are 
valid also for regions given by radial ranges ri < r < r2 not containing a symmetry 
center, so that S^ ' need not vanish at the extremes ri, r2. In particular, we can even 
apply the results proven in the lemmas to LTB models that lack a symmetry center. 
In these cases (all of which are elliptic) , the integral definition for quasi-local scalars 
([5]) may not converge, hence it is more convenient to define and compute niq, kg, 7iq 



first, and then use (14aH(14c| to obtain the local scalars m, k, H. 

As we show further below, a change of sign of kq at some r = rj, introduces specific 
constraints on 5'^'^^ that can affect the criterion for the existence of TV's and profile 
inversions. The possible combinations of glueing hyperbolic and elliptic regions with a 
parabolic region, either containing a center or not, force the constrain kq = to hold 
in a given radial range, which is rather artificial, hence such combined constructions 
will not be examined here. 

10.3. Mixed elliptic/hyperbolic construction. 

Since an elliptic region enclosing a parabolic or hyperbolic region containing a center 
is incompatible with the Hellaby-Lake conditions [17], we consider the opposite 
construction made by an elliptic region containing a symmetry center surrounded 
by an expanding hyperbolic exterior. This interesting mixed configuration can be 
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k. 



constructed by prescribing a given niqi > together with 

> for < r < Tf), elhptic region 

= for r — rb, interface . (92) 

< for r > ri,, hyperbohc region 

where we assume that fc'j and rn'j are continuous (at least C^) at r = r;,. Since 
k'i{rb) < and fcgi(r(,) = 0, then 



qi 



as r ^ r\, we have 5, 



(fe) 



— cx) for r < rt, 
oo for r > ri,, 



elhptic region 
hyperbolic region 



(93) 



which implies that S^ > must hold in the hyperbolic region at least in the 
neighborhood of r = r^. From section 6 and Lemma 9, S^"^^ < and (5'-'°' < 

(k) . J 

must hold everywhere in the elliptic region. However, the diverging of 5- in (931 



does not allow us to test the Hellaby-Lake conditions (C.7) and (C.8) and t o eva luate 
^(™) and S^'^^ at r = Tb by setting kqi = (or equivalently Xi = 0) in (22), (45&) and 
(63 1. Instead, since r = rj, corresponds to kqi — 0, the form of these functions in the 



neighborhood of r = r^ follows form power series expansions on Xi ~ kqi/niqi = 

1 



nqc{t 



1 _ Ik. 
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Oix,), 



SO that cC, r, J(™) and (5^=) in (22), (C.5), {45b) and (63) take the forms 




(94 a) 
(946) 

(95 a) 

(956) 

(95c) 
(95 d) 



(fe) 



where all the terms 0(xi) involve products of the form kqiSl , which are finite even of 
61 ' diverges. Therefore, ct'^^^, T and 6^"^'> change smoothly as kq passes from positive 
to negative in the interface r = rh, and simply take their corresponding parabolic forms 
around the inte rface, so that regularity there is fulfilled if the parabolic Hellaby-Lake 
condition (C.6) holds. Also, (5^™^ < necessarily holds in the interface and (by 



continuit y) a t lea st in its vicinity inside the hyperbolic region. From Lemma 9, the 
limits in (93) and {95 d\, it is evident that that 5^'^^ < holds for all times in the elliptic 
region, and in particular S^'^' —>■ — oo holds in the limit r —>■ rb within this region, while 
in the hyperbolic region we have J^*''-' — >■ oo as r — >■ r^, and thus (5^'^' > must hold at 
least in the neighborhood of r = r^, all of which is consistent with k'Jt^rb) < and 
kq{t,ri,) = holding for all t. 

As a consequence of the restrictions on (5'^'^^ described above, there cannot be 
density or spatial curvature profile inversions in the elliptic region, in which niq and 
kq have clump profiles (and thus m and k by Lemma 3). However, profile inversions 
of these scalars are possible in the hyperbolic region depending on the possibility to 
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attain the right combinations of behavior of (5^ and 5^ that yield these inversions 
(as displayed in Table 1). 

For r > rii but close to ri, we must have ^ > and i5 < 0, which corresponds 



to case (iii) in Table 1, so that (49 1 does not hold. For these values of r (close to 
rb) there cannot be a density profile inversion, but a spatial curvature inversion must 
have happened for some t < ti. Since fc^ < holds for all r > rb, then we must have 
kg — >■ — oo as i — >■ tbb (or L — > 0) for every r > rb, which means that S'^^^ < must 
have occurred for r close to rb in hypersurfaces ■^T[t] close to ttb, but S^ > for 
t — ti, hence the profile inversion at t < ti. 

For r > rb (but sufficiently far from rb) in the hyperbolic region we can have the four 
combinations of profile evolutions outlined by Table 1: 

• Case (iii). If S^ > and (5^ " < holds for all r we have the same situation 
as in the layers close to r;,: 5^™^ < and (5''^) > hold for t > ti and thus 
there are no profile inversions of ruq and kq for t > ti. In this case S^'^^ > for 
kg < implies that kg < holds for for all r, with kg ^ possibly as r — >• cx), so 
that density decreases (clump) and spatial curvature becomes more negative as r 
increases (void profile for negative kg). 

• Case (i). If 6^ (ri) — with ri > r^ and becomes negative for all r > ri, with 
Jj remaining negative for all r, then a density TV and density profile inversion 



will occur for some t > ti in the range r > ri if ( 49 1 holds 



Case (ii) . As above but ( 49 ) does not hold in the range r > ri , hence there is no 
density profile inversion for t > t^, but there is a profile inversion of kg. 

Case (iv). As above, but (5^ (r2) = with r2 > ri and becomes positive for all 
r > r2. In this case mgi has already a void profile {m' ^ > 0) for r > r2, and since 



( 49 1 holds this profile is kept for the whole evolution (though a density profile 



inversion must have occurred for some t < ti). 
Evidently, the initial value functions rrigi, kgi, 5 and (5- must also comply with 



the regularity conditions (C.7), which do not follow from the qualitative restrictions 
described above. 

11. Summary and conclusion. 

We have conducted a comprehensive and rigorous examination of the radial profiles 
of the main covariant scalars (density, spatial curvature and expansion, respectively 
denoted by m, k, H, see definition ([t])) for regular LTB models in full generality. 
The time evolution of these profiles has been considered, and in particular, we have 
addressed the issue of profile inversions in which an initial clump can evolve into a 
void (or vice versa). The necessary background material is the formalism of quasi- 
local variables within an initial value formulation (section 3), rephrasing in terms of 
this formalism the analytic solutions (Appendix B) and the Hellaby-Lake conditions 
(Appendix C) , and the relation between the radial coordinate and proper radial length 
(Appendix D). The quasi-local density (jUg), spatial curvature {kg) and expansion 
scalar {'Hq) are given by ^-(10), while the fluctuations 5'^^'^ follow from (13). 
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11.1. Summary of results. 

We have used thoughout the article the term "turning value" of any scalar A ( "TV of 
A" ) to denote a value r = rtv such that A' vanishes a,t r — rtv at a given hypersurface 
^T[t] marked by constant t (this definition also applies to the Ag). In general, the 
value rtv is different for each '^T[t], but it can be the same for all t (as for example 
a regular TV of R). Throughout sections 3-8 we have considered only LTB models 
or regions containing a symmetry center and with an open topology (i?' > holds 
everywhere: no TV of R), leaving the case with closed topology (regular TV of R) and 
special configurations for sections 9 and 10. The main results contained in 10 lemmas 
that were proven in these sections are summarized below: 

• Section 3. The clump/void character of radial profiles of local and quasi-local 
scalars, A and Ag, in radial domains containing a center was given in terms of 
the existence of turning values, TV's. Assuming that R' > holds, we proved 
the following lemmas: 

— Lemma 1: The signs of the radial gradient A' is related to the sign of the 
fiuctuation S^'^\ 

— Lemma 2: The necessary and sufficient condition for a TV of Ag at some 
r — rtv is given by S^^> = at r = rtv. In general, rtv — rtv[t] is different at 
different times (see figure 1). 

— Lemma 3: The existence of a TV of Ag in a radial domain is a sufficient 
condition for the existence of a TV of A in the same domain. 

These lemmas allow us to examine the monotonicity of scalars by looking at signs 
and the zeros of their fluctuations. As a consequence of Lemma 3, it is sufficient 
to examine the profiles of the scalars Ag (which satisfy simpler scaling laws) to 
know the profiles of the scalars A. 

• Section 4. Character of initial density and spatial curvature profiles compatible 
with the Hellaby-Lake conditions. Initial density must have a clump profile in 
parabolic and elliptic models or regions, while initial clump and void profiles 
are possible in hyperbolic models or regions. Initial curvature must have a clump 
profile in elliptic models or regions, but hyperbolic models or regions admit initial 
clumps and voids. Explicit conditions are given for each case. 

• Section 5. Formal definition of profile inversion and its relation with TV's 
of scalars. Profile inversions can occur with and without TV's. The following 
lemmas were proven: 

— Lemma 4: A profile inversion of Ag implies a profile inversion of A. 

— Lemma 5: Sufficient conditions for the existence of a profile inversion of Ag 
follow from the existence of TVs of Ag plus some extra requirements on the 
fluctuations S^^' . This type of profile inversion is illustrated by figure [ij 

— Lemma 6: Sufficient conditions for the existence of a profile inversion of Ag 
when there are no TV's of Ag. This type of profile inversion is illustrated by 
figure [2j 

Notice that Lemma 4 guarantees that the results of Lemmas 5 and 6 apply also 
to the local scalars A. The conditions in Lemmas 5 and 6 are independent of the 
Hellaby-Lake conditions, hence the latter place extra constraints when applying 
these lemmas to specific configurations. 

• Section 6. Density radial profiles and profile inversions are examined in detail. 
The following lemmas were proven: 
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— Lemma 7: Necessary and sufficient condition are given for tire 'clump to 
void' density profile inversion with density TV's in regular hyperbolic models 
or regions. Void to clump inversion is not possible. 

— Lemma 8: Necessary and sufficient condition for the 'clump to void' 
inversion without TV's in regular hyperbolic models or regions. Void to 
clump inversion is not possible. 

— Lemma 9: Density profile inversions, with or without density TV's, cannot 
occur in regular open elliptic models or regions containing a center (no TV's 
ofi?). 

The main result of this section (and possible of the whole article) follows as a 
consequence of these lemmas: 

the only density profile inversion (with and without TV's) compatible 
with absence of shell- crossings (the Hellaby-Lake conditions) is the 
'clump to void' inversion in hyperbolic models or regions. 



The specific conditions outlined by Lemmas 7 and 8 (equations (49) and (51)), 
which must be satisfied by hyperbolic models and region, are simple restrictions 
on their initial value functions (density, spatial curvature and their fluctuations) . 
It is evident that these are not weird or outlandish conditions, but reasonable and 
easy to prescribe. 

Section 7. Spatial curvature profiles and their inversions. The only inversion 
compatible with the Hellaby-Lake conditions is the 'clump to void' inversion 
in hyperbolic models or regions. These inversions only occur in the particular 
parameter cases when there is no inversion of the density profile. 

Table 1 provides the combination of initial conditions that allow for profile 
inversions of density and spatial curvature in hyperbolic models or regions. 

Section 8. Radial profiles and profile inversions of the expansion scalar. 
Hyperbolic models and regions allow for a regular evolution with and without 
profile inversions (conditions are provided for each case). However, profile 
inversions of the expansion scalar necessarily occur in elliptic models or regions 
containing a center. This is illustrated by figure |3] 

Section 9. Closed elliptic models. A TV of R occurs under regular conditions 
at the same value r = rtv for all times. We proved the following lemma: 

— Lemma 10: The TV of i? at r = rtv implies that rtv is a common TV of the 
local and quasi-local density, spatial curvature and expansion scalar. The 
converse is not true, as we proved in Lemmas 1-5 and 7-8 that TV's of these 
scalars can occur when there is no TV of R [R' > holds everywere). 

As an important corollary of this lemma, the fluctuations S^ ' do not change sign 
because of the TV of R. As a consequence, there are no profile inversions of 
density and spatial curvature, though profile inversions of the expansion scalar 
must occur as in open elliptic models. 

Section 10. Special configurations: 

— Regular hyperbolic and elliptic models with a simultaneous big bang. Pro- 
file inversions of density and spatial curvature (with or without TV's) cannot 
occur. Hyperbolic models allow for a void profile for all the time evolution. 

— Regions not containing symmetry centers. The results of all proven lemmas 
hold, the only caveat being the fact that (in general) we have S'-^^ 7^ at 
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the boundaries of the radial range (as opposed to 6^'^'> — strictly holding 
at the symmetry center). 

— LTB models with "mixed" kinematics . These configurations are made by 
glueing or matching combinations of parabolic, hyperbolic and elliptic regions 
(see [SH E2 EH SI])- The results of all proven lemmas hold. 

~ Elliptic region surrounded by a hyperbolic exterior. From Lemma 9 and 
assuming that Hellaby-Lake conditions hold, the elliptic region must have 
clump profiles for the density and spatial curvature, without inversions. 
However, void profiles and 'clump to void' inversions are possible in the 
hyperbolic exterior in agreement with Lemma 7 and the combination of initial 
conditions specified in Table 1. 

11.2. Relaxation of regularity conditions. 

In obtaining the results summarized above we have assumed that the Hellaby-Lake 
conditions [lOlEIllSH] (see Appendix C) hold for all the evolution times of the models 
(shell crossings are completely absent). Evidently, if we relax these conditions by 
demanding that they hold only for all t > ti for a given t = ti, then some of the 
results of the Lemmas are no longer binding. The most important example is Lemma 
9, which forbids density void profiles and 'clump to void' density profile inversion 
in elliptic models. This result follows from the fact that a density void profile is 



incompatible with the joint fulfillment of two of the Hellaby-Lake conditions (C 
t'bb — ^-nd t'oii ^ Oi where itt and icou denote the locus of the initial and collapsing 
singularity. However, if we demand i'^,j > to hold but not t^^ < 0, then a shell 
crossing singularity necessarily emerges for t k, t^b for all r, but not for t k, t,.^ii. 
Depending on the free parameters ruqi^ kqi, it may be possible to construct an elliptic 
model that is free from shell crossings for all t > ti for some ti > ibb [M], and that 
admits a 'clump to void' inversion and density void profiles in this time range. Since 
a dust source does not provide a good description of the physical conditions near 
an initial singularity in a cosmological model, it may be sufficient for a physically 
reasonable evolution to demand absence of shell crossings for all t > ti, provided one 
can justify that the dust description beaks down for t < ti. 

The relaxation of the Hellaby-Lake conditions, as described above, can also work 
for hyperbolic models or for special configurations. It can also allow for a density 
'void to clump' profile inversion, and it can also be set up so that an evolution free 
from shell crossings occurs in time ranges restricted to an early stage 1^,^, < t < ti. 
This latter case can be useful for studying void formation in the context of toy models 
of structure formation ( "spherical collapse model" ) , in which the collapsing stage of 
an over-density with positive spatial curvature (elliptic) is omitted, as the late stage 
of the dynamics is dominated by the virialization process leading to the formation of 
stable structures (hence there is no collapse stage because the assumption of spherical 
dust is no longer valid) i48j. 

11.3. Final discussion. 

The most important results in this article are those obtained in section 6 (Lemmas 
7, 8 and 9), section 9 (Lemma 10 and its corollaries) and the elliptic/hyperbolic 
mixed configuration in section 10. These results correct and provide full generality to 
the work initiated by Mustapha and Hellaby [Mj, as these authors only furnished a 
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restricted proof (together with particular and numeric examples) of the existence of the 
'clump to void' density profile inversion in hyperbolic models, without attempting to 
prove it for elliptic or parabolic models. They also claimed (mistakenly) that the 'void 
to clump' inversion was possible without shell crossings. Since they only examined 
density profiles, the results of sections 7 and 8 also extend and generalize their work 
to spatial curvature and the expansion scalar. 

The fact that the existence of density void profiles (with total absence of shell 
crossings) is only possible in hyperbolic models and regions is consistent with the 
conclusions of [SIl [311 123] • However, as opposed to these authors, we have obtained 
the analytic conditions for these profiles and also for their inversions in full generality, 
and without resorting to artificial configurations made by matchings regions in which 
the free parameters take special forms and radial gradients and density are likely 
discontinuous. Since the results of Lemma 9 are analytic and general, the elliptic 
models with void profiles examined in [37] must necessarily have shell crossings (which 
was not proven otherwise by the authors). Also, some of the configurations in the first 
and second column of Table 1 of [33] are listed as free from shell crossings. This is 
mistaken, as the Hellaby-Lake conditions are incompatible with elliptic and parabolic 
regions that contain a symmetry center and have density void profiles. 

Radial profiles of scalars in closed elliptic models are strongly constrained by the 
TV of R that occurs at a fixed value r — rtv, which (as proven by Lemma 10) is a 
common TV for the density, spatial curvature, expansion scalar and all initial value 
functions (but not the initial fluctuations). Since this TV implies a change of sign 
of the gradients m'^ and m' at r — rtv, it introduces relative under-dense or over- 
dense regions in thick spherical shells of dust layers, as opposed to a compact spherical 
region around the center. These shells can be considered as some sort of "topological" 
density clumps or voids (the same remark applies to spatial curvature or the expansion 
scalar), and are a product of the spherical S'^ topology of the rest frames ^T[t], and 
as such are not related to a profile inversion produced by the evolution of the models 
given generic initial conditions (as is the case for the TV's that were examined in 
sections 3-8). As a consequence, models of topological clumps/void profiles [^ I5D] 
are rather artificial constructions. 

It is important to emphasize that the radial profile of the expansion scalar, as 
examined in section 8, is much less restrictive to the existence of TV's and profile 
inversions than the profiles of density and spatial curvature. This fact is crucial 
to assess the existence and magnitude of a "back-reaction" term and an "effective" 
acceleration in the context of Buchert's scalar averaging applied to LTB models 
[2niini]i which is a relevant issue in current research involving LTB models. Also, as 
pointed out in recent literature [49] , spatial gradients of the expansion scalar (related 
to gradients of the binding energy) are the main generators of back-reaction (see 
also |20[ I21| ). and as such are important in providing a theoretical interpretation of 
cosmological observations that avoid introducing a dark energy source. 

Regarding the elliptic region surrounded by a hyperbolic external region, the 
fact that the latter may admit a density void profile (and a 'clump to void' 
inversion), without violating regularity conditions (Hellaby-Lake), is indeed a nice 
and unexpected result. This possibility leads to simple models in which a high 
density region undergoing collapse is surrounded by a large spherical shell in an 
intermediary transition scale, whose density is lower than that of the cosmic large 
scale asymptotic background. Configurations of this type could be compatible with 
the notion of "finite infinity" in which cosmic bound structures are approximately 
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regarded as asymptotically flat at such intermediate scales between the cosmic voids 
and the far cosmic background gHl ISO] ■ 

Finally, it is necessary to remark that we have conducted a theoretical study of 
radial profiles and profile inversions based on general analytic expressiond (avoiding 
the use of special cases based on narrow parameter specializations) . Hence, we have not 
taken into consideration the question of whether these radial profiles are compatible 
or not with the constraints that follow from actual cosmological and astrophysical 
observations. However, introducing specific profiles meeting these observational 
criteria is an important task that certainly requires a follow up article based on a 
full numeric approach, preferably using the Omega and Hubble parameters discussed 
in Appendix E. This article is presently under elaboration. 

Appendix A. Analytic solutions and the Hellaby Lake conditions in the 
conventional variables. 

The solutions of the Friedman-like field equation (pi) for each kinematic class take the 
following well known parametric form: 

Parabolic models or regions: E — 0. 

c(i-U-^^', i^-(2M)l/3,?^ (A.l) 

Hyperbolic models or regions: _E > 0. 

i? = — (cosh 7] - 1) , c{t- i, J = -^^ (sinh 77 - 77) , (A.2) 

Elliptic models or regions: E < 0. 

i? = 7— (1 - cos 77) , c{t - ibb) = n^pT^ iv - sin 7/) , (A.3) 

where itb = ^bb(?') is called "big bang time", as it marks the coordinate locus of 
the central expanding curvature singularity: R{t, r) = for r > 0. This free function 
emerges as an "integration constant" in the integration of (pi) . Notice that the locus of 
central curvature singularity is distinct from that of the center of symmetry R(t, 0) = 
(see Appendix Al of [38 J. Besides the kinematic class, LTB models that admit 
(at least) one symmetry center can be classified as "open" or "closed", respectively 
corresponding to the hypersurfaces of constant t being topologically equivalent to M.^ 
or S'^ (see Appendix D and Appendix A3 of [38j). 

Appendix B. The analytic solutions in terms of an initial value approach. 

The Friedman-like equation (p|) is equivalent to (10), hence the analytic solutions of 



the former (given by (A.1|-(A.3)) are equivalent to those of the latter, but given in 
terms of L and the initial value functions {ruqi, kqi}, which are related to the free 
parameters M and E by (11). Notice that {niqi, kqi} is an irreducible set of free 



functions, since Ri can be specified as a radial coordinate gauge. As we show below, 
the "bang time", tbb(''), follows as a function of 771^^, kqi. See Appendix E for an 
alternative set of initial value functions that can be related to the Hubble and Omega 
factors of a FLRW cosmology. 
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We express M and R in terms of ruqi and L by inserting (|6|) and ( 11 ) in ( A.l ). After 
re-arranging terms we get the following closed analytic expression for L 

2/3 

(B.l) 



L = 



1 + -^/2m~c{t-ti) 



where we are only considering expanding configurations {L increases for t > ti). The 
bang time follows by considering that Li — 1 and setting L = and i = tbb in (B.l ) 

2 



Ctbb — Cti 



3^27 



Cti 



^qi 



S'Hqi 



(B.2) 



Appendix B.2. Hyperbolic models or regions: kqi < 



We obtain the following implicit solution of the form t = t{R) by eliminating the 
parameter rj from the equation for R in (A. 2 1 and substituting in the equation for t: 

^3/2 

-j^c(t-t^,^)=Zh{R), (B.3) 

where R — {E/M)R and Zh is the function 

u h-> Zh{u) = u^/2 (2 + uf^ - arccosh(l + u). (B.4) 

We express then M, E and R in (B.3 1 in terms of ruqi^ kqi and L from (IgI) and (11). 
The result is 



yi c[t - ti) = ZhixiL) - Zh{xi), 



where 



Vi 



\kq^\'^' 



(B.5) 



(B.6) 



^qt 



f'qi 



Setting L^ = 1, t ~ ty^^ and L = in (B.5) and using (B.4t yields the bang time as a 
function of viqi and \kqi\\ 



c^bb — cti 



Vi 



(B.7) 



Appendix B.3. Elliptic models or regions: kqi > 



In this case, the implicit solution t — t[R) follows from eliminating rj from the first 
equation in (A.3) and substituting in the second one. The resulting implicit solution 
has two branches, an "expanding" one (0 < 77 < tt with i? > 0) and a "collapsing" one 
(tt < 77 < 27r with 7? < 0): 

\E?'^ 

-c{t-t^^) -- 



M 



Z^(R) expanding phase 

27r — Ze (R) collapsing phase 



where R = {\E\/M)R and Ze is given by 

u t-^ Zeiu) = arccos(l ~ u) ~ v}!"^ (2 - uf^"^ . 
Proceeding as in the hyperbolic case, we transform (B.8I into 

Ze(XjL) 



yic(t-ti) + Z^{xi) 



27r - Z^(x^L) 



expanding phase 
collapsing phase 



(B.8) 
(B.9) 

(B.IO) 
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where 



Vz 



^qt 



l3/2 



''qi 



(B.ll) 



It follows from (|B.10|)-(|B.9|) that L is restricted by < L < i„ax, where the maximal 

^,qi, characterized by L = and Hq = 0. 



2mqi/k, 

and L = in the expanding phase of (iBJOl) yields a bang 

,1 and L 



expansion is L,„ax — S/x^ 
Setting Li = 1, t — t 
time function, while t = t^^„ and L = in the collapsing phase yields the "crunch" 
time associated with the collapsing singularity. The maximal expansion time follows 
by substituting t = i^^x and L 



in either branch of (B.lOl. These times are 



given by 



cih 



ct„ 



ct. 



Cti 



Cty, 



Ze\Xi) 



= ci,. 



27r 
Vi 



cti 



= cti 



Ze(Xi 



27r — Z^ (xi 



(B.12) 
(B.13) 

("•^^) 
tma^ir) and i^„„ = icoii(?'), like t-t,i,{r), are not 
simultaneous (do not coincide with a '^T[t] hypersurface). For every comoving observer 
r = const., the time evolution is contained in the range ibb('') <t< t^^n^r). 

The parametrization of initial conditions given here in terms of {rriqi, kqi} is 
not unique. We could have used {mqi, Hqi} as initial value functions, so that 



Notice that (in general) t„ 



^ruqi — T-LZi follows from ( 10 ). An alternative description (which can be suitable 



for applications) is provided in Appendix E in terms of initial value functions that 
resemble observational parameters (the Hubble and Omega parameters of a FLRW 
cosmology) . 



Appendix C. Avoidance of shell crossings: the Hellaby Lake conditions. 

The restrictions on the free parameters {M, E, cibb} that guarantee an evolution free 
from shell crossings (thus complying with (24)) are the well known Hellaby-Lake (nec- 
essary and sufficient) conditions [40l Ell I46j given for each kinematic class by 



Parabolic and hyperbolic models or regions: 



R' >Q ^ 
Elliptic models or regions: 

± i?' > -^ 



{Af'>0, E'>0, C<0}. 



± M' > 0, ±C < 0, 



± 



2 E 



2ttM 



(C.l) 



>0, 



(C.2) 



where only expanding configurations are considered in (C.l) and the ± sign in (C.2) 



accounts for the fact that R' < occurs in elliptic models whose that admit a second 
symmetry center (see section 9 and Appendix A3 of [SB])- The equal sign holds only 
at symmetry centers and at values of r where R' = 0. 
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In order to express the Hellaby-Lake conditions in terms of the quasi-local scalars 
in the initial value formulation M we examine the sign condition of the analytic forms 



of r in (21) and (22) f35j. A comparison with the forms (C.I) and (C.2) follows by 



expressing the gradients Af, £" and ctj^j, in terms of our initial value functions by 
means of 



M' _ 

ctL 



Ri 



1 



r(™) 



SR^Ri 
iR[/Ri 



c{ti - ibb) 5. 



^(™) _ ^(k) 



(m) 



Ey_ 
'e 

2 6^ 



Ri 



(m) 



3^2^ 



qi 



n 



- c{t, - tbb) S, 



qi 



:(»") 



Ak) 



parabolic 



2 ' 



hyperbolic and elliptic 



(C.3) 



(C.4) 



(C.5) 



where we differentiated both sides of (11 ) 



(13) specialized to t = ti, while itt in (C.5 



B.2) and (B.7) with respect to r and used 



) is given by (B.7) with Zh{xi) or Ze{xi 



respectively, for hyperbolic and elliptic models. Considering (C.3) and (|C.4[)-(C.5), 



the Hellaby-Lake conditions for each kinematic class are: 
• Parabolic models or regions. 

- 1 < ^f"' < 0, 



(C.6) 



which from (13) implies m'^ < 0. It is evident from (C.3) and (C.4) that (^C.6) 



"qi 



-Lara 



implies tL/R'i < and M' > 0, and so it is equivalent to the Hellaby-Lake 
conditions (C.l ). 



• Hyperbolic models or regions. 



> 



jtI™) 



> -1. 



(C.7) 



wher e tL is given in terms of our initial value functions by (C.5). Notice, from 
(186), that all regular hyperbolic models comply with S'-'^^ > —2/3 for all their 



the Hellaby-Lake conditions (C.l). 
Elliptic models or regions. 



time evolution. It is evident from (C.3) that (C.7) are completely equivalent to 



CtL 



iR'^/Ri 



<0, 



CtL 



with 



ct' 



- U, 



-(m) 



3R^/Ri 



3^(^-) 
2 * 



>0, 



c(™) 



c(icoii - ^bb) + 



>-l, 



CtL 



(C.8) 



(C.9) 



where we us ed c(^ eoii - ti) = c(teoii - ibb) - c{ti - tbb), with itt and teou given by 
(B.12) and ( B.14[ ), and the equal sign holds only at a symmetry center. Notice 
that two of the three conditions in (C.8) are sign conditions on the gradients of 



the coordinate locus of the central singularity {tL, ^Lu)- This fact has not been, 
apparently, noticed in previous literature [H [HI EH US] . 

(t We examine the Hellaby-Lake conditions in the parametrization {rriqi, kqi\, as it is straightforward 
to re— phrase all expressions in terms of the initial conditions {'Hqi, Oi} defined in Appendix E 
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Conditions (C.8) imply the following necessary (but not sufficient) condition for 

m 



s-(m) 



^^f ^ > 0, 



(C.fO) 



which, from ( 18c), implies that (5^™^ — (3/2)6^'''' > necessarily holds for all times 
if (24 1 holds. Since condition ( C.10[ ) implies 

l + 3(<5f"'-5f^)>0,, (C.ll) 



this (also necessary but not sufficient) condition follows from (C.8) too. It is 



straightforward to show from (13), (B.12), (B.14), (C.3), (C.5| and (C.9) that 



conditions (C.8) are equivalent to the Hellaby-Lake conditions (|C.2| 



Appendix D. Proper radial length and a radial coordinate gauge. 

As mentioned in section 3, the radial coordinate r is not intrinsically covariant. We 
provide in this Appendix the conditions guaranteeing that the dependence of scalars 
on r is qualitatively analogous to their dependence on the proper radial length. We 
also discuss the choice of radial coordinate, given the existence of a radial coordinate 
gauge freedom in LTB models by virtue of the invariance of the LTB metrics (fl]) and 
(20 1 with respect to an arbitrary rescaling r = r{f). 

Radial radial rays are spacelike geodesies of the LTB metric [39l 



whose affine 

parameter is the proper radial length, which (if there is a symmetry center r = 0) is 
given by the function: £[i] : M+ — > M such that 



mir) 



. Ax 



ViTe 



dx 



LT 



yr 



f^qi ^i 



: dx. 



(D.l) 



where we used (11) to eliminate E in terms of kqi. Since £ must be a non-negative 



and monotonously increasing continuous function, the regularity condition (24) must 
hold together with the following extra regularity conditions 



,.R^ 



>o, 



1 






>o, 



so that a zero of R'^ (if it exists) must be a common same order zero of -y/l — kqiRf. 
Since the zeroes of R[ and R' are common "HS", "3^ , this condition automatically follows 
if (24) holds. The converse is not true, as (D.2) can be violated even if (24) holds. In 



this latter case, the radial coordinate would be ill defined (as well as the proper radial 
length), and also a surface layer singularity would arise at the common zero of R[ and 
1-kg^R^ gH [461147]. 

For regular "open" models whose ^T[t] are topologically equivalent to M.^, or 
for any LTB region in which R' > and 1 — kqiRf > hold everywhere without 
violating ( 24 ) and ( |D.2 ) , Ri can be prescribed as any monotonously increasing function 
complying with Ri{0) — and R'^ > for all r. Evidently, the simplest choice in these 
cases is 



Ri = Rq r, 



(D.3) 



where i?o is an arbitrary constant characteristic length scale. The gauge (D.3) is a 
popular choice in the literature [TOl HSl [El [17j [18] , not only due to its simplicity, but 
because by setting r = R^/ Rq radial dependence becomes dependence on a fiducial 
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value of an invariant quantity that has a clear physical and geometric meaning (R). 
Also, the choice of Rq provides a physical length scale for the radial coordinate. 

In "closed" models whose '^T[t] are homeomorphic to S'^ (see section 9) the 
function Ri cannot be monotonously increasing (and so (D.3) cannot be used). In this 
case there are two symmetry centers, so that Ri{0) 
value" r = rtv must exist so that i?',(rtv 
r = rtv must be also a zero of 1 



0, hence a "turning 



fvqiR 



R'{ct,rtv) = 0, but because of (D.2| 
^ of the same order. Regular closed models 



cannot be parabolic or hyperbolic because kgi < holds for these models, hence (D.2 1 
would be violated [HJ \^ [37] . As a consequence, all regular closed models must be 
elliptic (and thus collapsing to a second curvature singularity). Since R'^ changes sign 



range rtv < r < r^ (see [Ml [Ml [4TJ]T 



at rtv, the regularity condition (D.2 1 implies that -y/l — kqiRf < must hold in the 
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The fact that Tiq behaves in ( 10 1 as a Hubble scalar of a FLRW dust universe suggests 
the definition of a sort of the following "Omega" quotient 



n 



n- 1 



2mq 

kq 
I Lq 



2m„ 



2to„ 



ZiTfti-i rCfi 



2m„ 



kqiL 



kqiJ-J 



2m„ 



2ra„ 



kqiL 



(E.l) 
(E.2) 



Lq n,q ^noqi 

This quantity has been used in various articles [13 [T71 [3S] and is only one of the 
possible generalizations of the Omega parameter for LTB models (see [M] for a proper 



discussion on this issue). The definitions (E.1|-(E.2) lead to 
2m„, 



VL, 



2mn 



'^li 



2m„ 



n,. 



1 



Krt 



SO that the following scaling laws hold: 



'^li 



2m„ 



n,~{n,-\)L 



n- 1 



Kf] 



{n. 



1)L 



n, - {n, -i)L 



riq — riqi 



i3/2 



(E.3) 



(E.4) 



(E.5) 



Notice that prescribing as initial condition a given sign for fi^ — 1 completely determines 
the kinematic class (parabolic, elliptic or hyperbolic) for all times (just as prescribing 



^qi) 



From (E.2) and (E.4) we have 



If f7i = 1 => r2 = 1, parabolic (E.6) 

IfO<0, <1 => 0<r2<l, hyperboHc (E.7) 

\i VLi>l ^ Cl>l, elliptic (E.8) 

Also, it is important to remark that irrespective of the kinematic class we have for 
every LTB model: f7 — > 1 as L — > (near a curvature singularity) . For all hyperbolic 
models or regions fii is bounded for all choices of rriqi, kqi, hence: f2 — >■ as L — )■ oo, 
whereas for elliptic models or regions f2 — )> cxd as i — > L^^^. 

We can always eliminate rriqi, kqi by means of (E.3) in terms of the initial value 
functions Cli, T-Lqi. Since the latter are roughly equivalent to Omega and Hubble 
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factors, or at least they could tend at a given asymptotic limit to these parameters in 
a FLRW background (see [35j for a discussion on this point), it could be more intuitive 
to prescribe fi^, Hqi as initial conditions than prescribing ruqi, kqi as in section 4. Since 
the parabolic case is trivial [Vti = 1, 'H^j = 2mqi), we provide below the hyperbolic 
and elliptic analytic solutions (B.5I and (B.IO) rewritten in terms of Cti, Hqi'. 

• Hyperbolic models or regions: (li — I < 0. 



W 



W,= 



c{t-t,) = 

1/2 



W-Wi 

T~Lai 



VLi + {l-VLi)L 



il/2 



n,. 



i-a 



1 



; arccosh 

2(1-^03/2 \n 



2(1 - ^,Y/^ 

2 



arccosh ( — — h 1 — 2L 



• Elliptic models or regions: £7^ — 1 > 
W -W^ 



1 



HqiC{t - ti) 



TrCliini - 1)-^/'^ - w - w. 



fli f2L 

W = —^ arccos U^ + 1 - 2i - 

2(r!,- 1)3/2 \^^ 

Wj — — ^ arccos \ ^ 1 

2(17,-1)3/2 \^^ 



expanding phase 
collapsing phase 

1/2 



n, - [a, - i)L 
a-i 



Li/2 



1 



a-1 

Setting Li — \^ t — t^,^ and L = in these expressiond yields the bang time 

itqi 

The time t = t^^^ and t — t^^„ in elliptic models, respectively associated with maximal 
expansion {L = i,„ax — ^i/i^i ~ 1)) and the collapsing singularity (L = in the 
collapsing phase) are 



(E.9) 

(E.IO) 
(E.ll) 

(E.12) 

(E.13) 
(E.14) 

(E.15) 



C = ct 



ttO," 



2nqd^,~ir/^' 



t^coll Ct^i-, + 



ttQj 



nq,[n,-i]y^ 



(E.16) 



The relation between Qi, Hqi and the conventional free parameters is 



a; = 



2M 



2M + E R, 



'Hqi — 



[2M + ER, 



,1/2 



R 



3/2 



(E.17) 



while using (13 1 and (E.3) we can express (5^ and S^ in terms of fli, Hqi as 



sr(™) _ 



3i?' 



n' 2H' 



qi 



n, 



Hn 



Ak) _ 



Ri 



Ci,-i 



2H\, 

Hn 



qi 



(E.18) 



where in (E.17) and (E.18) the function Ri can be specified as a choice of radial 
coordinate gauge. 
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Appendix F. Local scalar representation. 

Besides p given by (l3|, otlier covariant objects associated witli LTB models are the 
expansion scalar, 8, the Ricci scalar of the space slices, ^TZ, plus the shear and electric 
Weyl tensors, aab, Eat 

^ ~ „ 2R k . 2{ER)' ,^^^ 

(Tab - ^iaUb) - (e/3)/i„b = S S'^^ E'^'' - u^uaC^^-" = E E-\ (F.2) 

where Va = /i^Vb and C'""^ is the Weyl tensor, while E"'' = h"'' - 3rfr]\ with 
Tj°' = ^/h/^S!^ being the unit tangent vector along the radial rays (orthogonal to w" 



and to the orbits of S0(3)). The scalars £ and E in (F.2) are 
R R' 






R R' 



^-E^^^''---,P+§- (F-3) 



Since LTB models (as all spherically symmetric spacetimes) are LRS (locally 
rotationally symmetric [42 ), they can be completely characterized by covariant scalars. 
Considering pi, (|F.1|, (F.2) and (F.3), a choice of scalar representation follows from 



the local "fluid flow" scalars [37l [38] 

{p, e, 37^, S, f }, (F.4) 

whose evolution equations completely determine the dynamics of LTB models in 
the fluid flow or "1+3" approach [l^ [331 [33], and thus provide an alternative 
approach to that based on the analytic solutions of (pi). Our usage of quasi-local 
scalars provides an alternative complete representation of covariant scalars given by 
{nig, Hg, kg, (5("), 5(«), ^W} [371 li. 
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